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FOREWORD 


The development of the analytical and test programs presented herein 
was sponsored by the National Aeronautics and Space Administration, Lewis 
Research Center, under Contract NAS3-21708. The NASA Project Manager was 
Mr . L . J . Kiraly . 

Principal United Technologies Research Center (UTRC) participants in 
the contract activity were Dr. A. V. Srinivasan, Mr. D. G. Cutts, 

Dr. S. Sridhar and Mr. J. Zucker. Dr. Srinivasan was the Principal Inves- 
tigator and Program Manager with primary overall responsibility for the 
development and application of the analyses and coordinating the test phases 
of the program. Dr. Sridhar provided the bulk of the support in the analy- 
tical aspects of the program and was primarily responsible in formulating 
the mathematical models and obtaining numerical solution through a series of 
computer programs he developed. Mr. Cutts provided the bulk of the support 
in the test programs and was primarily responsible for outlining the tests, 
designing the fixtures and instrumentation, directing the tests and data 
reduction. Mr. Cutts was assisted by Mr. John Zucker throughout the test 
programs . 

This report is the final documentation of (a) analytical modeling of 
nonaerodynamic sources of damping in engine blades, and (b) testing per- 
formed on blades and blade— like components to measure the levels of damping 
available in these components. 


i 



LIST OF FIGURES 


Page 


Blades and Other Test Components Used in the Program . . .107 

Twisted Plate Mounted on Centrifuge Arm 108 

Twisted Plate Installation for Centrifugal Test 109 

Turbine Blade/Platform Damper Vibration Test Set Up . . .110 

Hysteresis Loop for Material Damping Ill 

Typical Blade Element . . m 

Stress Distribution Across Element Thickness Ill 

Material Damping Constant Measurement Test Hardware . . .112 
Assembly of Centrifugal Test Rig for the Twisted Plate . .113 

Twisted Plate Mounted in Bench Test Fixture 114 

NASTRAN Finite Element Model for the Twisted Plate .... 115 
Block Diagram of the Driving and Recording Systems 
Used to Obtain Vibration Decay Data for the Twisted 

Plate on the Centrifuge 116 

Variation of Modal Damping (Loss Factor) with 
Rotational Speed for Titanium Twisted Plate in 

Partial Vacuum (5 torr) 117 

Decay Characteristics of Twisted Plate Vibrating in 

Its First Three Modes at Zero Speed in Vacuum 118 

Decay Characteristics of Twisted Plate Vibrating in 

Its First Mode at Speed in Vacuum (Two Samples) 119 

Decay Characteristics of Twisted Plate Vibrating in 


Its Second and Third Modes at Speed in Vacuum (5 torr) . .120 
Instrumented Fan Blade Showing "Root" Strain Gage . . . .121 
Instrumented Fan Blade Showing A.S.M.T. Strain Gage . . .122 
Test Set-Up For Fan Blade Material Damping Investigation 


at Elevated Temperatures 123 

Variation of Measured Modal Damping with Air Pressure 

for Fan Blade 124 

Variation of Measured First Mode Frequency Ratio with 

Air Pressure for Fan Blade 125 

Material Damping in a Typical Fan Blade 126 

Variation of Loss Factor with Root Stress and Soak 

Temperature 

Stress Amplitude Decay Curves for a Fan Blade Vibrating 
in Air In Its Second Bending Mode for Two Temperature 

Conditions 128 

Stress Amplitude Decay Curves for a Fan Blade 
Vibrating in Air In Its First Torsion Mode for 

Two Temperature Conditions 129 

Typical Composite Material Test Piece 130 



LIST OF FIGURES (Cont’d) 


Figure Title Page 


26 Composite Material Cantilevered Beam Vibration 

Test Set Up 131 

27 Composite Material Damping 132 

28 Spanwise Distributions of Outer Fiber Major Principal 

Stresses for Twisted Plate 133 

29 Twisted Plate Forced Response in Partial Vacuum (10 Torr) . 134 

30a Motions of Blade with Slip at Root 135 

30b Analytical Model of Blade with Slip at Root 135 

31 Test Fixture and Root Loading Device for Fan Blade 

Root Damping Tests 136 

32 Fan Blade/Fixture Test Assembly for Root Damping 

Investigation 137 

33 Circuit and Switching Arrangement for Strain 

Gage Load Bolts 138 

34 Decays of Blade Response with an Applied Root 

Loading of 576N 139 

35 Decays of Blade Response with an Applied Root 

Loading of 2879N 140 

36 Damping at Fan Blade Root 141 

37 Variation of Root Vibratory Stress with Root Loading . , . 142 

38 Variation of Blade Response Frequency with Root 

Normal Load for Various Excitation Levels 143 

39a, b Schematic of Shroud Rubbing 144 

40 Analytical Model for Shroud Rubbing 144 

41 Shroud Loading Device for Investigation of Fan Blade 

Shroud Damping Mechanism 145 

42 Shroud Loading Mechanism 146 

43 Block Diagram of Basic Excitation and Recording Systems 
Used to Obtain Strain and Acceleration Response Data 

From Test Items Excited by the Electro-dynamic Shaker . . 147 

44 Test Set-Up for Shroud Damping Investigation 148 

45 Fan Blade Response Characteristics from a 2g Level 

Exploratory Sweep - Shrouds Loaded with 102N 149 

46 Fan Blade Response Characteristics from a 2g Level 

Exploratory Sweep - Shrouds Loaded with 150N 150 

47 First Flap Mode Response to a Slow Frequency Sweep at 

Constant Input Acceleration Levels 151 

48 Damping Due to Rubbing at Shroud for a Typical Fan 

Blade - Variation of Response Decay with Initial Stress. . 152 

49 Damping Due to Rubbing at Shroud for a Typical Fan 

Blade - Variation of Response Decay with Shroud Load and 
Initial Stress 153 

50 Shroud Damping Effects on Blade Response (First Flap 

Mode ^-280 Hz) 154 


iii 



LIST OF FIGURES (Cont’d) 


Figure Title Page 

51 Blade Response in the Above-Shroud Torsion Mode 155 

52 Wear Areas on Shroud Load Platens 156 

53 Friction Test Assembly (Upper Arm Raised) 157 

54 Friction Test Assembly 158 

55a Typical Signals from Acceleration and Force Transducers . .159 

55b Typical Measured Friction Force - Slip Loops I 59 

56 Variation of Equivalent Friction Coefficient During 

Sinusoidal Motion With Normal Load, Frequency and 

Maximum Relative Velocity 160 

57a Schematic of Turbine Blade and Platform Damper 161 

57b Analytical Model for Platform Damping 161 

58a Damper Force vs Blade Displacement at Platform 162 

58b Blade and Damper Spring Displacements Over One 

Cycle of Vibration 162 

59 Turbine Blade Test Piece Geometry 163 

60 Damper Stress Response Signals for Varying Input 

Acceleration Levels 164 

61 Baseline Conditions for Turbine Blade 164 

62 Platform Damping Test Results (Large Damper) 165 

63 Platform Damping Test Results (Small Damper) 166 

64 Turbine Blade Platform Damping Characteristics 

(Small Damper) 167 

65 Influence of Damper Stiffness on Platform Damping 168 


iv 



LIST OF TABLES 


Table Title Page 


I Titanium Cantilever Beam Test Results 96 

II Summary of Twisted Plate Test Results 97 

III Summary of Data for Fan Blade Material Damping Tests, ... 98 

IV Summary of Root Damping Test Data 99 

V Summary of Shrouded Blade Response Test Data 100 

VI Summary of Friction Test Results 101 

VII Summary of Turbine Blade Response Test Results 102 

VIII Summary of Turbine Blade Platform Damping Test Results 

("Small Damper") 103 

IX Summary of Turbine Blade Platform Damping Test Results 

("Large Damper") 104 

B-1 Test Components and Other Hardware Designed and Fabricated 

for the Test Program 105 

B-2 Equipment Used in Test Program 106 


V 



NOMENCLATURE 


area of beam cross-section; surface area of element 

amplitude, nondimens ional, of forced response 

measure of slip amplitude for platform damper 

constants, in Section 3.2 

linear viscous damping coefficient 

energy dissipated during one cycle of vibration 

as a subscript denotes damper location for platform damper 

modulus of elasticity of material 

frictional force, nondimensional, at an interface 

as a subscript, denotes friction 

slip amplitude, nondimensional 

acceleration due to gravity 

as a subscript, denotes hysteresis 

thickness of uniform cantilever beam 

mass moment of inertia of fan blade about root 

cross-sectional area moment of inertia of beam 

as a subscript, denotes element number 

damping constant, for material damping 

modal stiffness; also used as cantilever beam frequency parameter 
in Section 3.1 

length of beam, and fan blade 

as a subscript, denotes lower surface for material damping 
calculations 



NOMENCLATURE (Cont’d) 


M moments; bending, frictional, and vibratory 

m modal mass 

N normal load at an interface 

n damping constant for material damping calculations; also used 

as subscript to denotes mode number in Section 3.2 

P amplitude of external excitation 

p amplitude, nondimensional , of external excitation 

Q ratio of response amplitude to external excitation amplitude, 

"magnification factor" 

q constants in Section 3.2 

i 

R effective radius to the rubbing surface for a dovetail root 

t time 

U strain energy stored during one cycle of vibration 

X forced response; also tip deflection of blade 

a sine and cosine components of the amplitude magnification factor 

3 stiffness ratio for shroud damping and mass ratio for root damping 

Y,5,6 phase angles 

e small dimensionless parameter; also used to denote strain in 

Section 3.2 

5 linear viscous damping constant, non-dimensional 

K correction factor for dovetail geometry 

)i coefficient of dry friction 

p mass density 
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T time parameter, nondimens lonal, for platform damping 

<p modeshape; also used as phase angle 

0) circular frequency 
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SUMMARY 


The potentials of various sources of nonaerodynamic damping in engine 
blading are evaluated through a combination of advanced analysis and testing. 
The sources studied include material hysteresis, dry friction at shroud and 
root-disk interfaces as well as at platform type external dampers. A limited 
series of tests were conducted to evaluate damping capacities of composite 
materials (B/AL, B/AL/Ti) and thermal barrier coatings. Further, certain 
basic experiments were performed on titanium specimens to establish the 
characteristics of sliding friction and to determine material damping con- 
stants J and n. 

The test components used in this program included part span shrouded fan 
blades, turbine blades, a twisted plate and uniform beams. Fixtures, loading 
devices, excitation and instrumentation systems were designed for use in each 
series of tests taking special care to highlight the mechanism under study 
and to minimize fixture participation. Measurement methods adopted for use in 
this program varied over a wide range requiring modifications to conventional 
procedures. 

The series of tests indicated that contributions to damping from 
material hysteresis is negligible for the titanium alloy used in this program. 
Further, the effect of root damping diminishes rapidly to negligible levels 
at projected operating speeds of a typical advanced fan with a dovetail type 
of root structure. A simulated thermal environment typical of a fan blade 
application showed essentially no changes in material damping characteristics. 
Similarly, the series of tests in which a twisted titanium plate was spun 
indicated essentially no changes in material damping levels. 

Damping at shroud interfaces is most likely to be the only nonaerodynamic 
source available at all speeds to limit vibration of shrouded fan blades. The 
nature of this damping has been found to be a complicated mechanism involving 
both micro- and macroslip conditions. The stress amplitude at slip resonance 
has been found to be proportional to shroud load and excitation levels whereas 
off resonance it is proportional to excitation levels only. The boundary con- 
ditions at the interfaces appear to change during vibration resulting in a 
very narrow region of contact. 

The mechanism of damping due to external dampers such as platform type 
dampers is, like shroud damping, inherently nonlinear. Depending upon the 
loading on the platform, the response characteristics change from a nearly 
sinusoidal signal to a stick-slip signal to a response emitting audible noise. 
Contributions to stress attenuation can come from damping due to rubbing as 
well as due to a change in mode shape brought about by the extent of stiffness 
of the damper. Test data showed changes in mode shapes occurring as the 
principal test parameters (excitation level, normal load) varied. 
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Thermal barrier coatings appear to have no influence on damping of 
turbine blade vibrations. Tests on composite material specimens indicate the 
problem of repeatability and suggest a broader program of tests involving a 
larger number of specimens and tests. 

A unique fixture design has been developed from which fundamental sliding 
friction characteristics of rubbing surfaces can be measured. The laboratory 
set-up is such that friction-displacement loops can be obtained from which 
work done/cycle of motion as well as coefficient of friction can be measured. 
For the titanium pieces with tungsten carbide surface treatment used in this 
program, the coefficient of friction was found to be between .3 and .5 depend- 
ing upon the normal load and independent of velocity. 

Although all the tests were conducted on single blades and appropriate 
mathematical models were developed to represent the several mechanisms of 
damping, procedures to apply this data to predict damping levels in an assembly 
of blades have been developed and discussed. 
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1 . INTRODUCTION 


Vibration induced fatigue failure of blades is of continuing concern to 
the designer of aircraft engines. The emphasis on improved engine performance 
under the necessary constraints of minimum weight and satisfactory life 
requires that vibration levels be kept low. Further, certain important de- 
sign considerations require a thorough understanding of the structural 
dynamic characteristics of blades. These design considerations include (1) 
blade life prediction methods which use vibration amplitudes in their calcu- 
lations, (2) allowable frequency margins which need to be justified on the 
basis of the intensity of resonant stresses computed at low integral orders, 
and (3) accurate prediction of susceptibility to aeroelastic instabilities. 

The analytical techniques used in flutter and resonant stress prediction rely 
on accurate determination of modal characteristics of blades. These charac- 
teristics include not only the natural frequencies and mode shapes but also 
damping available in each mode. Contributions to damping in a vibrating 
blade arise from aerodynamic sources as well as nonaerodynamic sources. The 
latter include contributions from material damping, friction damping due to 
rubbing at interrupted interfaces (shroud to shroud, root to disk), and fric- 
tion damping due to an external device such as a turbine blade platform-type 
damper. 

Damping, in the context of this report, refers to energy dissipation 
capacities of the component undergoing vibration. A measure of energy dissi- 
pated per vibrational cycle is commonly used to report the extent of damping 
present in a system. In this report, the extent of damping is reported as 
the loss factor, which is defined as the ratio of energy lost per cycle of 
vibration to maximum vibrational energy stored in the appropriate mode 
divided by 2'it. Where such a measurement is inadequate or impractical, the 
extent of damping present is measured indirectly and reported as attenuation 
of strain in the vibrating components due to the influence of the damping 
mechanism. 

Material or hysteretic damping refers to energy dissipation due to many 
complex mechanisms within a material, when a volume of the material is sub- 
jected to cyclic stresses. Thus, material damping is always present in a 
vibrating blade. Friction damping at an interface, however, depends on load- 
ing, roughness of surfaces, level of external excitation, slip amplitude, 
geometry of the contacting components, etc. For example, the root structure 
of a jet engine blade is typically of a dovetail, pin or firtree design, the 
extent of friction damping being different in each design. In the case of 
shrouded blades, untwisting of the blades under centrifugal loading brings 
neighboring blades into contact at shrouds, resulting in a complex joint 
at the interface. Determination of the precise nature of boundary conditions 
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at such an interface continues to be an unresolved problem defying both 
analysis and measurement. The condition at the interface may vary from a 
fully stuck to a freely slipping condition, the extent of friction damping 
depending on the nature of slipping motion. External dampers are frequently 
used to reduce stress levels in turbine blades. The two common types of 
dampers are the so-called blade-to-ground and blade-to-blade dampers. During 
engine operation, centrifugal forces bring the damper into contact with the 
blade platform, the extent of friction damping depending on the nature of 
slipping motion at the contacting surface. 

The primary objective of this combined analytical and experimental pro- 
gram was to gain an understanding of the potentials of nonaerodynamic sources 
of damping in reducing vibration levels in jet engine blades. A critical 
survey was conducted of the available literature dealing with vibration damp- 
ing, in general, and damping in jet engine blades, in particular. Though 
literature on the subject of damping is extensive, relatively few studies 
appear to be directly applicable to jet engine blades. The list of references 
given at the end of this report (Refs. 1 through 48) is not meant to be exhaus- 
tive and only those studies which are of direct relevance to this program are 
included in the list. Based on this survey, a series of analytical and experi- 
mental efforts were identified for detailed study. 

In this program, single blades and blade-like components were analyzed 
and tested to obtain the contribution from each mechanism of damping. The 
components included a part-span-shrouded fan blade, a high pressure turbine 
blade, a blade-like titanium twisted plate, uniform cantilever beams of ti- 
tanium, and strips of composite material. In addition, small titanium speci- 
mens were used in rub tests aimed at determining the nature of dry friction 
characteristics at an interface. 

No single analytical or experimental approach appeared feasible for 
evaluation of all sources of damping, due to the diversity of mechanisms in- 
volved. Analytical efforts encompassed the mathematical modeling and analysis 
of each mechanism, leading to the solution of many independent problems. 
Experiments were conducted with the primary goal of highlighting the particu- 
lar mechanism under study and every care was exercised to minimize contribu- 
tions from other sources. The environments that needed simulation included 
centrifugal loading and, to a limited extent, thermal condition. Special fix- 
tures were designed and fabricated for each series of tests. In the case of 
the twisted plate centrifugal tests were conducted in an evacuated centrifugal 
test facility. 

This report summarizes (1) the principal features of mathematical modeling 
of the several damping mechanisms, (2) the results obtained from a series of 
tests leading to the measurement of damping associated with each mechanism, (3) 
a comparison between results from analyses and corresponding tests, and (4) 
implications of the results of this study to blade design against resonance 
and flutter. 
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2 . TECHNICAL APPROACH 


The fan blade and the turbine blade that were chosen for use in this 
program are shown in Fig. 1 along with other test components which were used 
to obtain additional supporting data. The survey of literature (Refs. 1 
through 48) indicated that various analytical and experimental approaches 
were available for the determination of damping in such components. These 
approaches differed considerably in the level of sophistication, ease of appli- 
cation, and probable time needed to adapt and implement for the component 
under study. Since no single approach appeared uniquely suitable and advan- 
tageous, decisions to consider more than one approach were inevitably made. 

The final choices in each case depended on later developments, trends observed 
in tests, and scheduled time and cost. Brief discussions of the technical 
approaches adopted during the program are provided below. 


2.1 Analytical Approaches 

The mechanisms for which analytical approaches were sought can be broadly 
divided into two classes; (1) damping caused by material hysteresis, and (2) 
damping caused by friction forces induced at rubbing surfaces. 

2.1.1 Material Damping 

Material or hysteretic damping is a material property which is measured 
by the energy dissipated during cyclic strain in the material. Accurate measure- 
ments show that cyclic stress-strain always has a hysteresis loop (Refs. 1 and 
2) . The analytical approaches for estimating material damping can be broadly 
divided into (1) formulations which represent the energy loss per unit volume 
per cycle of vibration in terms of stress amplitude and two material damping 
constants (Refs. 1 and 2), (2) formulations which model the hysteresis loop 
by a nonlinear stress-strain law for use in a forced response analysis (Refs. 

5, 6, and 7), (3) formulations based on a nonexponential law of decay (Ref. 

8) , and (4) formulations which express the stress-strain law in terms of a 
complex modulus; storage and loss modulus (Ref. 10). The first two formula- 
tions appeared to be suitable for application in the present program and were 
selected for detailed study. 

Thus, in one approach, the hysteresis loop is modeled by a nonlinear 
stress-strain law. This law is introduced into the equations of motion of a 
vibrating component and a forced response solution is obtained for comparison 
with experimental results. Apart from the inherent difficulties of nonlinear 
analysis, such an approach is not readily applicable to components having com- 
plex geometries. The alternate approach starts with a representation of the 
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energy dissipated at a point in the material in terms of two material damping 
constants and the vibratory stress amplitude at the point. Total energy dissi- 
pated is obtained by summing the energy dissipated at all points in the com- 
ponent. The procedure calls for a knowledge of vibratory modal stress dis- 
tribution, which can be determined accurately even for components such as jet 
engine blades, having complex geometries. In view of its generality, the 
latter approach was adopted in this program. Subsequently, the approach using 
the nonlinear stress-strain law was used to develop a procedure for determining 
the two damping constants. 


The ratio of total energy dissipated (D) to the maximum strain energy, 

(U) , during one cycle of vibration, is used as the measure of material damping 
for a component. Thus, the component "loss factor" p is defined as 




JL ^ 

2ti TI 


( 2 . 1 ) 


2.1.2 Friction Damping 


Depending on conditions at an interface between two contacting components, 
energy dissipation due to friction may involve different mechanisms. The 
analytical approaches for estimating friction damping can be broadly classi- 
fied as follows: (1) macroslip approach (Refs. 13 through 17), and (2) micro- 
slip approach (Refs. 23 through 25). In the macroslip approach the entire 
interface is assumed to be either slipping or stuck. The friction mechanism 
is either replaced by an equivalent linear viscous model or assumed to be 
governed by some form of Coulomb's law of dry friction. The effect of fric- 
tion damping is obtained by determining the forced response of the component. 
The analysis involved is relatively straightforward and the justification for 
widespread use of this approach is its effectiveness at predicting the actual 
response. However, there is some question as to the validity of the macroslip 
approach when the interface has a large area, and is subjected to a nonuniform 
load distribution. In the microslip approach a relatively detailed analysis 
of the stress distribution at the interface is carried out, typically via a 
finite-element procedure. The extent of local slip, not necessarily throughout 
the interface, between pairs of contacting points is determined by applying 
Coulomb's law of friction to the normal and tangential stresses. A detailed 
knowledge of interface slip dynamics can be obtained from this approach. 
However, the computational effort needed to modify, iteratively, the inter- 
face stress distribution subsequent to local slipping can make the microslip 
approach unacceptable in a practical application. 


In the context of jet engine blades, the macroslip approach was found to 
be suitable for modeling the friction mechanism and was therefore adopted in 
this program. Friction damping in jet engine blades arises from slip at the 
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root, shroud rubbing, and rubbing due. to a platform damper. The problem in 
each case can be described, in a broad sense, as one of determining the effect 
of friction forces, at some location along the span, on the vibratory motion 
of the blade. The mathematical model in each case is derived by using a modal 
approach. The governing equation of motion for a blade can be written as 
(Ref. 46) 


[M]{q} + [c]{q} + [K]{q} = {f}. (2.2) 


where M, K represent the mass and stiffness matrices. The vector F includes 
external excitation as well as any friction damping forces. All other damp- 
ing forces, such as material damping, are included in the equivalent viscous 
damping matrix C. The free vibration solution of Eq. (2.2) provides natural 
frequencies, mode shapes, and modal masses. By defining a modal matrix 

w = >. . . . n. 

assuming that, 

{q} = [l>]{x(t)}, 


(2.3a) 


(2.3b) 


and using the orthogonality property of the mode shapes, Eq. (2.2) can be 
transformed into 




where modal damping has been assumed. 

Assuming that external excitation is harmonic, the main interest is in the 
case where the excitation frequency is near one of the blade natural fre- 
quencies; that is, the near resonance condition. In fact, it is precisely 
at this near resonance condition that damping plays a crucial role in limiting 
the otherwise high blade stresses. Since the thrust of this program is 
toward evaluating the influence of damping on blade responses, attention is 
focused on a single chosen mode of blade vibration. Thus, the matrix Eq. (2.4), 
can be reduced to 

mx + cx + kx = P COS ojt + Fjj, (2.5) 
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where m, c, and k are modal quantities associated with the mode under 
consideration. P is the excitation amplitude, to is the excitation fre- 
quency, and accounts for friction damping forces. Equation (2.5), with 
appropriate modifications and extensions, forms the basis of study for each 
of the friction damping mechanisms, considered. The precise form of F^j, of 
course, depends on the particular source of friction damping (i.e., root 
slipping, shroud rubbing, and friction due to platform damper). 

An assumption which is implicit in the use of Eq. (2.5) is that the 
^^Tction damping does not drastically alter blade resonant frequency and 
mode shape, so that attention can be focused on a single chosen mode of 
vibration. The validity of this assumption is questionable if the analytical 
model is to cover the entire range of conditions at an interface, the two 
extremes being the fully stuck (no slip) condition and the freely slipping 
(zero normal load) condition. However, if the vibratory environment is such 
that the condition at an interface remains close to one of the extremes, then - 
the assumption made above is valid, and use of Eq. (2.5) is justified. 

In this report, the measure of friction damping is given as the increase 
of loss factor and/or attenuation of blade response over a "baseline” or 
reference condition. The loss factor is defined in Eq. (2.1) (i.e., propor- 
tional to the ratio of energy dissipated to maximum strain energy, during one 
cycle of vibration) . The baseline condition is simply a condition where the 
vibratory characteristics of the blade (tip amplitude, blade stress) are known 
For example, the condition of no slip at the interface (i.e., no friction 
damping) can serve as a baseline condition. 

The principal features of the computer programs developed during the 
course of this investigation are presented in Appendix A. 


2.2 Experimental Procedures 

The choice of the experimental procedure to measure the contribution of 
a specific mechanism to the overall damping in a jet engine blade is governed 
by a number of requirements derived from considerations of the design and 
operational aspects. The requirements are: 

a. the test item is generally small with a complex shape; 

b. the test item is cantilevered but with the possibility of relative motion 
at the root; 

c. damping to be measured in a single mode; 
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d. some very low damping levels to be expected; 

e. enviromental considerations to include a centrifugal force field, 

elevated temperatures; 

f. rubbing friction effects at root and shroud to be studied; 

g. air damping effects to be eliminated or accounted for. 

This list is not exhaustive but represents the major considerations used 
when planning the various aspects of the experimental procedures, i.e., (1) 
fixture design, (2) method of excitation, (3) environmental simulation, and 
(4) measurement method. 

2.2.1 Fixture Design 

In tests involving the measurement of damping, enormous care needs to be 
exercised in fixture design so that fixture participation, both in resonant 
response and damping, can be minimized. In the present program, special 
fixtures were designed and fabricated with the goals of providing as realis- 
tic an environment as possible and highlighting the particular damping 
mechanism under study. 

The basic fixture arrangement used for measuring material damping in 
blade-like components was a massive root block set in a rigid (within the 
test frequency range) vacuum chamber and mounted to a slip plate so that 
vibration could be applied in the horizontal direction. This method was 
chosen because it provided a convenient way of applying adequate excitation 
force to the component in an evacuated environment and also for providing 
means for loading the shrouds and dovetail bearing surfaces of the fan blade. 
The block, vacuum chamber, slip plate, and shaker armature formed a massive 
root restraint for the test component. 

The component used in the centrifugal tests was a titanium twisted plate 
having an integral massive root bayonet fitting which locked into the end 
of the rotating arm; see Figs. 2 and 3. A special locknut arrangement was 
devised to ensure that the loading in the bayonet fitting did not vary with 
rotational speed. 

The fixture for platform damping investigation was basically a heavy 
clamp to hold the root of the turbine blade and a mechanical means provided 
to allow a damper to rub on the platform of the blade; see Fig. 4. The 
assembly was mounted on top of an electrodynamic shaker. An evacuated 
environment was considered unnecessary because of the expected small blade 
amplitude and high damping levels from the external damper. 
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2.2.2 Methods of Excitation 


The external excitation on a blade or blade-like component can be provided 
by a hammer impulse, an initial tip deflection, a sinusoidally varying tip 
force from a noncontacting magnetic driver or a sinusoidally varying root 
displacement. For each specific test, the method of excitation and the 
design of the fixture had to be chosen on the basis of the force levels re- 
quired, the ease of operation and control, and applicability to the measure- 
ment method to be employed. In this program all excitation (except that used 
in centrifugal tests), was sinusoidally varying root displacement provided 
by an electrodynamic shaker. In the centrifugal tests, the twisted plate 
was excited by a noncontacting electromagnetic driver located at the tip. 

2.2.3 Environment Simulation 


Since this program is concerned with evaluation of nonaerodynamic sources 
of damping, every effort was made to minimize, if not eliminate, air damping 
effects. Test components were enclosed in a specially designed chamber con- 
nected to a vacuum pump. A centrifugal test facility was used to apply cen- 
trifugal loading to a twisted plate in an evacuated condition. In an operat- 
ing engine, centrifugal forces lead to interfacial loading between contacting 
components. Since all friction damping tests in this program were bench tests, 
the interfacial loading was applied and controlled by specially designed 
loading mechanisms. Tests conducted to study the effect of temperature on 
damping of a fan blade called for the simulation of a typical operational 
thermal environment. This was done by using blowdrier heater elements to 
provide radiation heating. 

2.2.4 Measurement Methods 

Measurement methods fall into three general categories, (1) decay measure- 
ment, (2) frequency response analysis, and (3) direct measurement. The 
decay measurement involves the initial excitation of the test item and an 
analysis of the subsequent response signature when the excitation force is 
suddenly removed. The form of the transient is a decaying sinusoid for 
response in a single mode. If the envelope of this response is considered, 
the loss factor can be determined from a plot of log amplitude versus time 
(log decrement - Refs. 33 and 34). A variant of this method is to average 
over a large number of cycles of the quasi-transient behavior of a system 
subjected to random excitation (random decrement) (see Ref. 34). 

The difficulty in decay measurements using the log decrement method is in 
producing a transient with a single frequency content. The problem is minimized 
if an initial sinusoidal excitation of a single mode is employed and the modes 
are well spaced (Ref. 7). More general transients, however, can be analyzed by 
using fast Fourier transform techniques to obtain frequency response functions. 

If this response (both magnitude and phase) is examined close to a modal frequency. 
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the modal damping can be obtained by a ''Nyquist” curve fit around the resonance 
The procedure for performing this has been automated in computing systems avail 
able on the market (see Refs. 35-37, 38). 

The frequency response method involves driving the test item through its 
resonance either by a sinusoidal sweep or dwells at discrete frequencies. 
Response amplification as well as phase angle between output motion and input 
force or motion are examined (Ref. 20). One measure of loss factor for single 
degree of freedom systems is the value of the inverse of transmissibility 
(Ref. 39). Another, of course, is by way of measuring the bandwidth at 
half power points (Ref. 40). 

The frequency phase method described in Ref. 41 enables the loss factor 
to be determined by measurement of the input frequency as well as the phase 
between response and exciting force. These measurements are made at a fre- 
quency slightly off resonance. The calculation is independent of the gain 
setting of the equipment and hence the problem of maintaining constant input 
level is obviated. 

The direct measurement method (Refs. 42 and 43) requires the system to be 
driven sinusoidally at constant amplitude at or near resonance and the energy 
input is measured. Since measurement is made at constant amplitude, the 
energy dissipated per cycle is equal to the average energy input per cycle. 

The loss factor can be found from the ratio of this value and an independently 
determined value of the maximum kinetic energy in a cycle. The estimate of 
the latter may be made analytically. 

The Co and Quad (see Refs. 44 and 45) method is similar to the modal 
analysis (FFT) method in that it separates the response vector (in transfer 
function form) into its real and imaginary parts and uses a Nyquist plot to 
determine the modal damping. However, the input in this case is a sine sweep 
at constant level. Alternatively, the damping factor may also be obtained 
from the frequencies of the maximum and minimum amplitudes of the quadrature 
component about a resonance (see Refs. 44 and 45) 

Each of the above methods has its own advantages and disadvantages rela- 
tive to a particular application. One or more of the above methods were used 
in making measurements of damping in this program. In certain cases, certain 
modifications were made to these methods and will be described in the appro- 
priate sections of this report. 
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3. MATERIAL DAMPING 


The structural configuration of a jet engine blade varies from a nearly 
uniform, thick, cambered blade of low aspect ratio with little or no twist to a 
thin, long, twisted and bent aerofoil whose camber, thickness, and cross-sectional 
area may vary nonuniformly from root to tip. A complex blade geometry subjected 
to gas and centrifugal loading leads to a complex stress distribution during vi- 
bration. For an accurate analytical estimate of material damping (which is a fun- 
ction of vibratory amplitude), vibratory stress calculation must necessarily be 
made taking full account of blade geometry. An experimental approach may pose 
its own problems, such as the need to eliminate fixture and air damping effects and 
the inherent difficulties of measuring small amounts of damping. 

Studies dealing with material damping can be found in Refs. 1 through 12. 
Presented in this chapter, are the results of analysis and testing of a titanium 
cantilever beam. In addition results from tests conducted on composite material 
specimens and coated turbine blades are discussed. 


3.1 Determination of Material Damping in a Vibrating Blade 


Due to the complex geometry of the blade and the resulting complicated stress 
distribution during vibration, it was found that an effective approach for estim- 
ating material damping is the one where energy dissipation is represented in terms 
of stress amplitude and damping constants; see remarks in Section 2.1,1. An out- 
line of the analytical procedure is given below. 

For a vibrating blade, the elastic response at a point in the material is 
characterized by a hysteresis loop as shown in Fig. 5. The area enclosed by the 
loop is a measure of the energy dissipated per unit volume per cycle of vibration 
and is called specific damping energy. The specific damping energy is assumed to 
be of the form 

D = Jcr'^ i J = (3.1) 

where a is the maximum principal stress, o is the fatigue strength (Ref. 2), J 
and n are material constants. The energy dissipated in the i element (in the 
sense of an "element" in a finite element program) is given by 

Dj = J f cr^dVj = JAj f cr^dz, (3.2) 

V| z 

where v^, A^, and h. are the volume, area of the upper/lower surface, and thick- 
ness, respectively, of the element (see Fig. 6) . 
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(3.3) 


The strain energy in the element is given by 


^ 2E ^ 2E 


where E is the Young’s modulus of elasticity 

Because of the complicated geometry of a typical fan blade, three cases may 
arise, depending on the stress distribution across the thickness of the element, 
as shown in Fig. 6b. 

In Case 1, the stress is constant across the thickness. The other two cases 
have linear stress distributions with the subscripts u and £ denoting the upper 
and lower surfaces, respectively, of the element. For the three cases, it can 
be shown that Eqs. (3.2) and (3.3) yield the following: 


Case 1 


0i = JAjhiO=.", 


•Ji ■ 2E 


' AjhjO-j 


(3 . 4a) 


Cases 2 and 3 




TTtr 


r_ n+i _ n+i 

°'u_ 

Ob+C7i 


' 6E 


-3.-3 

Ob ± cr_£ 


(3.4b) 


where the minus signs apply to Case 2 and the plus signs apply to Case 3, and 
overbars denote absolute values. 

Here it is assumed that modal stress distribution of the vibrating blade is 
available via a finite element analysis. The material damping loss factor for 

the blade can be written as X Dj 

I J 

- 27T Jy. (3.5) 

i 

where the summations are carried out over all elements in the finite element model 
of the blade. 

A computer program, MATDMP, was written to compute the loss factor for a 
vibrating blade. The program accepts as input the maximum principal stress on 
the upper and lower surface of each element, the area and thickness of the element, 
the material constants, E, a , J*, n, and magnitude of a reference stress. It is 
noted that the loss factor for a vibrating blade is associated with a specific 
reference condition; for example, the maximum stress at the root is to be, say, 

60 MPa. The program calculates D. and U. for each element using the appropriate 
formulas in Eq. (3.4) and uses Eq^ (3.5) to calculate the loss factor. 
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It is noted that for the special case of n=2, the loss factor for the blade is 
independent of the stress distributions and reduces to the simple expression 



(3.6) 


The procedure given above for calculating material damping in a vibrating blade 
is quite general and can be used for components having complex geometries. It 
is noted that modal stress distribution is a standard output from finite element 
programs such as NASTRAN. However, the accuracy with which the loss factor is 
determined depends critically on the accuracy of the material damping constants 
used. The approach developed to obtain the damping constants is presented below. 


3.2 Determination of Damping Constants J and n 

Though the representation of material damping in terms of damping constants 
has been discussed in the literature (Refs. 1 and 2), reliable values of the 
constants appear to be unavailable. For example, the damping constants J and n 
tabulated in Ref. 2 for some materials are based on tests conducted at very low 
frequencies (0.33 Hz). There is some question as to whether such information 
obtained at essentially static test conditions can be utilized in a dynamic anal- 
ysis. Therefore, a need exists for a method of determining the material damping 
constants accurately over ranges of frequencies expected in application. 

The several methods proposed in Ref. 3 for determining the damping constants 
require the use of a specially designed test fixture. The methods are based on 
measuring one or more of the following quantities: energy input into the fixture; 

torsional deflection of a calibrated shaft within the fixture; and strain distrib- 
ution along a cantilever test specimen. It was felt that a major disadavntage of 
these methods is the need for a rather complex test setup with its inherent suscep- 
tibility to experimental errors. In the present work, a method of determining J 
and n from a simple test setup, involving vibration of a cantilever beam, is pre- 
sented. The required analysis is outlined below. 


The hysteresis loop, which characterizes material damping, is modeled by a 
nonlinear stress— strain law. The model is introduced into the differential eq- 
uation of motion of a vibrating uniform cantilever beam. The forced resonant 
response to a harmonic excitation is calculated. The solution for the response, 
which is a function of the excitation and damping, leads to a method of calculat- 
ing J and n from test results. 

The stress— strain law chosen is a modification of the one proposed in Ref. 5, 
and can be written as 

a- = cr^ ’ (3.7a) 
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(3.7b) 


°-N 




(3.7c) 


where o' is the nominal stress and is the contribution from material hysteresis. 
N H 

The upper signs refer to the ascending part of the loop, the lower signs to the 
descending part, ^ is the strain and is the strain amplitude. It can be shown 

that the area of hysteresis loop yields Eq. (3.1) (i.e., specific damping energy). 


The equation of motion of a uniform cantilever beam undergoing small-amplitude 
vibrations and subjected to harmonic base excitation is, 


__a^y __a y ^ '^h --2 -- 

El— + — ;- = paJ COS cut, 

ax'* ^ dx^ 


(3.8) 


where overbars denote dimensional quantities. Here y(x,t) is the deflection, p the 
mass density, A the area of cross section, I the cross-sectional area moment of 
inertia, p the base excitation amplitude, o) the excitation frequency, and the 
bending moment due to the presence of material hysteresis. The equation of motion 
can be nondimensionalized by introducing the following quantities: 


x = xL, y = yL, p = €pL, 



9 E” ' J /_h_\ 
+ \L/ 


+ [(yo"± y")^-2(yo'')*] . 


(3.9) 


where the cross section is assumed to be rectangular, h is the thickness, L the 
length, a prime denotes differentiation with respect to x, and y is the amplitude 
of vibration. The expressions for the parameter e, and the bending moment can 

be obtained by using Eq. (3.7) and the strain-displacement relationship of simple 
beam theory. The procedure is straightforward and the details are omitted here, 
for the sake of brevity. The final result is a nonlinear differential equation, 

3"^y . <3^y . 2 

T + — ^ + € = cpoj cos wt • (3.10) 

d>r dr ax^ 

As e 0, the linear free vibration problem is recovered. The solution of 
Eq. (3.10) is sought in the form 

y (x.t) = S <^„(x) ( 3 . 11 ) 
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Substituting 
leads to 


where cj) are linear eigenfunctions and u are yet to be determined. 

Eq. (3.5!i) into Eq, (3.10) and using the^othogonality property of cj) 

I 

“n “'on ‘Jfi = “ T [/ +C 2„pa)* cos (Ktl , (3.12a) 

In 0 

where 

Sn ’ °2n'/'^n • (3.12b) 

0 0 

and oa are the natural frequencies, 
on 

Due to the presence of the nonlinear hysteresis damping term, an exact 
solution of Eq. (3.12) cannot be obtained. However, by using the method of 
averaging (Ref. 47), a first-approximation to the solution can be obtained. As 
the interest is focussed on the resonant response in a chosen mode of vibration, 
the subscript n can be dropped and the solution written in the form 

u = Q (t) COS [cj^t + = a cos y. (3.13) 


The amplitude and phase are assumed to be functions of time, in contrast with a 
problem of linear vibrations where a and 0 are constants. 

The essence of the method of averaging is in the assumption that, if the 
nonlinearity is small (i.e., e is small), then the amplitude and phase may be 
treated as constants within a cycle of vibration . Application of the method 
transforms Eq. (3.13) to a set of first-order differential equations 


a = 


2C| Ojo 


[B| + Cg sin (y-wt) , 


® " Sc.La [*i • 


(3.14) 

(3.15) 


where to has been assumed to be near to , and A^, are the first Fourier 
coefficients in the expansion 

I 

f 4>" M dx = A (a) + I A (a) cos ny + X (a) sin n y. (3.16) 

Q n ' n 

By using the expression for M i.e., Eq. (3.9), Eq. (3.10) and Eq. (3.13), the 
Fourier coefficients can be shown to be, after some algebra, 
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(3.17) 


“■8 c ^ 

= 01 ' =s=N')"dx. 

0 

The problem of obtaining the first-approximation to the steady-state response 
reduces to the problem of solving two nonlinear algebraic equations 

a*^'' + p sin = 0 , (3.18) 

q € 0 - a + q a^"' + p cos \p = o , (3.19) 

3 2 


where 

'p = y-cji , 

CJ = (Jq ( 1 + € cr) , 

A *^2 2 

P = E7 % 


8 C 3 € 

2 Cjf 

... 2 


■ 7TC, ’ 

c. 

^3 “ 0 ■ 


From Eq, (3.18), it is seen that the maximum amplitude, §, for a given excitation, 
occurs when simjj = -1, so that. 


A n-i 

a 



(3.20) 


The linear eigenfunction of a cantilever beam is 


</> = c[(sin kx-sinh kx) + c (cos k x - cosh k x) 


(3.21a) 


where 


A 


C = 


COS k + coshk 
sin k - sinh k 


(3.21b) 


the characteristic number k is given by. 


cos k cosh k + I = 0 



(3.21c) 


and the normalizing coefficient C is chosen such that (f)(1) = 1. 
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In base excitation tests, it is convenient to measure the strain at some 

point on the beam, say at the root. The relationship between the root strain, 

£ , and the tip deflection, y, , can be shown to be 
o ■'L 

2 ^ 

'rr «0 ’ =-k c C . (3.22) 

n C., “ 

Since the objective is to develop a method of determining the damping 
constants from test results, an equation relating root strain to base excitation 
amplitude is required. By using relationships defined among the various quan- 
tities during the course of the analysis, it can be shown that Eq. (3.20) re- 
duces to 


where 


C 


5 


J 




(3.23a) 


12 C3 e"''l 

TTCg ' 2" (n + l) h 


(3.23b) 


Eq. (3.23a) can be rewritten as 

in = in (Cj j ) + (n-i) in €q . (3.24) 

Thus by plotting (p/L) vs e on a log-log scale, the damping constant n can be 
determined. Equation (3.23^ can then be used to determine J. 

The test data required for determination of the damping constants consists of 
excitation levels and corresponding responses for a uniform cantilever beam. A 
computer program DMPCON was written to generate a linear least-square fit to the 
test data, and calculate the damping constants.. The program accepts as input the 
geometry of the specimen (length, thickness), material properties (elastic modulus, 
mass density), test conditions (mode number, location of the strain gage on the 
specimen), number of test points, and a table of base excitation levels vs strain. 
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3.3 Material Damping Tests 


The series of tests that will be discussed in this section focussed on the 
determination of the characteristics of damping inherent in the titanium alloy 
(8-1-1). Thus, tests were performed (1) on a uniform cantilever beam to determine 
damping constants J and n, (2) on a twisted plate on the bench as well as in a 
centrifugal force field to determine damping, (3) on a fan blade to determine its 
material damping characteristics. All tests were conducted under evacuated con- 
ditions to minimize, if not eliminate, air damping effects. 

3.3.1 Uniform Cantilever Beam 

The objective of this test was to determine the material damping constants 
J and n for a typical titanium alloy in bending and at typical blade response 
frequencies. 

3 . 3 . 1 . 1 Test Approa ch 

In Section 3.2, a relationship between the maximum bending strain and peak 
input acceleration (sinusoidal forcing function) at the root of a uniform can- 
tilevered beam is derived in terms of the constants J and n. The log-log plot 
of the peak excitation amplitude versus the corresponding maximum root strain 
has been shown to be a straight line from which the constants can be determined. 
The approach chosen therefore was to tightly clamp a beam test piece at one end 
and subject it to harmonic excitation measuring the strain at a point close to 
the root for a range of peak excitation amplitudes and in each of its first 
three modes. The root strain was derived from the measured strain using the 
standard beam functions (Ref. 48). 

The test piece, as shown in Fig, 7a, was designed and built incorporating 
the following criteria some of which were in accordance with Ref. 39. 

a) Titanium 8-1-1 alloy to be used, 

b) Single cantilevered beam configuration, 

c) The clamped length of the beam to be at least one tenth the beam length 
and have a thickness of at least three times the beam thickness to min- 
imize clamping effects, 

d) The aspect ratio to be at least 6:1 to minimize Poisson's ratio effects, 

e) The root radius to be equal to the beam thickness to minimize the stress 
concentration factor, 

f) The uniformity of thickness and width along the beam length to be better 
than 1%, 

g) The first three resonant frequencies to be in the range 60 to 1200 Hz, 
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h) The piece to be machined from bar stock with successively reduced cuts 
to minimize residual stresses, 

i) Surfaces to be finely ground. 

A special clamp (see Fig. 7b) was designed to enable the test piece to be 
rigidly held in the vacuum chamber of the basic fixture assembly described in 
Appendix B. A torque of 27 m.N was found to be sufficient to clamp the test 
piece, as no change in beam response (vibrating in first mode in air) was ob- 
served as the torque was increased further. The input control accelerometer 
was mounted directly on the clamp near the root of the beam. 

3 . 3 . 1 . 2 Instrumenta tio n 

The instrumentation system used is basically the same as that discussed in 
detail in Appendix B with the exception that the X-’^ plotter was used only to 
record steady state strain levels (unfiltered, rms) with no horizontal sweep 
time base involved. The signal from the input control accelerometer (VSS Type 101, 
SN105) was amplified 20 dB prior to entering the servo control unit during exci- 
tation of the first beam bending mode so that the required low acceleration levels 
(g<.l) could be applied. The response of the blade was obtained from a strain 
gage (Micro-Measurements Type EA-15-062AK-120) mounted 17.2 mm from the root. The 
gage lead-in wires were routed and epoxied in such a way as to minimize any damp- 
ing or mass effects. 

3 . 3 . 1 . 3 Test Resu lts 

A fixture check out confirmed that there were no fixture resonances near 
the blade modal frequencies of interest. There was, however, a high frequency 
component noted in the input acceleration signal. This was eliminated by the 
use of a 50 Hz bandwidth tracking filter inserted in the excitation circuit 
between the signal conditioner and the servo control unit. 

Owing to the very low damping in the beam (p 'v .0002) the maximum response 
amplitude of the beam was difficult to obtain. The signal generator instability 
generally exceeded the half power bandwidth of the response. However, by care- 
ful tuning procedures using a Lissajous figure displayed on an oscilloscope and 
repeated attempts, a reliable value for the peak response was obtained. When 
the peak value of the response was considered to have been obtained, the strain 
gage level was read off the digital voltmeter and also recorded on the X-Y 
plotter. The excitation period was noted when the’ response ellipse axis on 
the scope was veritical. 

All tests were performed with the pressure inside the chamber reduced to 
15 torr. The final set of response data for each of the three beam modes is 
shown in Table I. These data were reduced using the expressions given in 
Appendix B and using the following equipment settings and other data: 
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Input accelerometer conditions setting 


41.8 mv/g 


Strain gage and lead in wire resistance 120.5 U 

Gage factor 1.995 

Strain gage translator setting 3.6 

Assumed modulus of elasticity 124 GPa 


3.3.2 Twisted Plate - Bench and Spin Tests 


Using the analysis of Section 3.1, the loss factor for a bladelike structure 
having a complex shape may be computed from the maximum principal stresses in 
a given mode. As a first step in applying this method to an actual fan blade, a 
uniform, twisted plate was chosen. The test objectives were to measure the forced 
response of the plate and to determine its loss factors both on the bench and 
with the plate subjected to a centrifugal force field. 

3 . 3 . 2 . 1 Te st Approach 

The analysis of Section 3,1 requires the determination of the principal 
stresses thorughout the component. However, if the responses are modal then the 
stress at one point can be used to scale the modal stress distributiom. Therefore 
the test method chosen was to harmonically excite the plate which was held rig- 
idly at its root in a massive block and measure the strain response using miniature 
strain gages at its root for various levels of forcing. The loss factors were 
estimated from root strain transient response curves resulting from abrupt cut 
off of the harmonic excitation, i.e,, in an essentially single mode response. 

The application of an overall steady force field was obtained by mounting the 
plate on the end of a centrifuge arm and spinning it. For this series of tests, 
the harmonic excitation was provided by means of an electromagnetic driver located 
at the tip of the plate. 

The test specimen was a titanium alloy 8-1-1 plate of aspect ratio, 3. 

The dimensions were 15.24 x 5.08 x .25 cm with the plate linearly twisted to 
give a twist of 30.5° (actual) relative to the root. Similar design and fabric- 
ation guidelines to those used for the uniform beam were used in this case (see 
Section 3.3.1). However, the root attachment of the centrifuge arm required a 
special design. The final design shown in Fig. 3, incorporated the following 
features : 

1. A root mass that was large compared with that of the plate. 

2. A minimum of mating surface area at the connection between the root and 
centrifuge arm. 
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3. Preloaded mating surfaces using a special locknut arrangement (torque 
to 680 Nm) such that the load remained constant over the range of test 
speeds , 

4. Integral plate and root bayonet fitting. 

A steel shim, 0.1 mm thick, was glued near the tip of the plate to permit 
magnetic excitation, see Fig. 2. The magnetic driver was attached to the arm by 
means of a steel tube surrounding the test piece. The test assembly in the cen- 
trifuge is shown in Fig. 8. The plate tip radius was 72.59 cm. 

For bench tests, the plate was mounted in the fixture/ vacuum chamber as 
shown in Fig. 9 and described in Appendix B, using a clamp adapter. This clamp 
was mounted on the massive steel root-block and held the bayonet fitting of the 
plate by means of four large bolts. The root chord of the plate was positioned 
vertically. 

The plate was instrumented with two longitudinal strain gages (Micromeasure- 
ments Type EA-06-062AK-120) mounted 0.64 cm from the root and edges of the plate 
as shown in Fig. 2. These positions corresponded with the root edge elements used 
in the finite element structural analysis model shown in Fig. 10. The length of 
wire on the actual plate was minimal with the wire routed along the radius and 
then up through a hole in the root. The wire was epoxied in place. 

3 . 3 . 2 . 2 Instru menta tion 

The instrumentation system used in the centrifugal tests is shown in a block 
diagram in Fig, 11. The plate was excited by means of a small electromagmetic 
driver located near the tip of the plate (see Fig. 2). Power for the driver was 
provided via slip rings by an MB Electronics power amplifier and the frequency 
was controlled by a B&K sweep oscillator. A Fluke counter provided an accurate 
readout of the input signal period. The plate strain gages were wired through 
slip rings to a pair of strain gage translators. These provided a calibrated 
signal to a narrow band pass filter and amplifier. The signal with high and low 
frequency noise removed was recorded on a B&K level recorder operating in the 
peak mode. Decays were obtained on this recorder and loss factors determined 
using the method described in Appendix B. The recording system used for the 
bench testing is essentially the same as the system described in Appendix B 
except that the narrow band filter was not required since the noise level was 
considerably lower. The excitation system used is described in Appendix B with 
the amplifier in the input circuit to enable the very low (< 0.02g) acceleration 
levels to be applied. The input control accelerometer (VSS type 101, S/N 105) 
was mounted directly onto the thickened root (bayonet fitting) of the plate test 
piece as shown in Fig. 8. For both tests, all output voltage rms readings were 
taken from a digital voltmeter and that level was noted on the charts. 
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3 . 3 . 2 . 3 Test ing and Res u lts 


Prior to testing in the centrifuge, the arm was carefully balanced using 
finely adjustable counter balance weights. Testing consisted of exciting the 
plate at each of its first three natural modes, namely first bending, second 
bending and first torsion modes, in air (at 0 rpm only) and in partial vacuum 
(5 torr) at rotational speeds up to 1508 rpm. The procedure used was to bring 
the plate to resonance, as indicated by a steady maximum response, and abruptly 
cut-off the excitation, thus obtaining a record of decaying signals from strain 
gages mounted near the root. The magnet had to be repositioned to obtain maximum 
response in each mode. The system used was found to be susceptible to noise 
from the slip rings and interference from neighboring rigs. However, by repeat- 
ing the tests at each test point a number of times, reasonably consistent results 
were obtained. In all, data were recorded from a total of 120 tests in a par- 
tially evacuated environment of 5 torr. A summary of the results obtained from 
these tests is given in Table II, and actual test data are plotted in Fig. 12. 
Sample decay curves are given for tests performed at zero speed in Fig. 13 and 
at 1000 rpm in Figs. 14 and 15. 

Bench testing of the plate consisted of exciting the plate in air and in 
partial vacuum in each of its first three modes with a base excitation, noting 
the peak response of the plate and then obtaining a decay curve after an abrupt 
cut-off of the excitation. The input levels needed to obtain adequate strain 
levels in the first mode were very low (< .02g) and required a signal amplifier 
in the input circuit. Even with this technique, the range for which values could 
be obtained was very limited. The second mode was difficult to tune because of a 
possible fixture interaction, but above a certain input level, consistent results 
were obtained. These results are summarized in Table II. Owing to small asymmetry 
of the plate, the first torsion mode was also excited but in this case the response 
to various inputs could not be realistically compared with analysis and so only 
the loss factor data are given. The values for loss factor in each of the modes 
from bench test are also plotted in Fig. 12 and provide a good comparison with 
values obtained from the centrifuge test. 

3.3.3 Typical Fan Blade in Ambient and Thermal Environments 

Using the analysis of Section 3.1, the loss factor for a component vibrating 
in a given resonant mode may be computed from the maximum principal stress dis- 
tribution. In order to verify the theory, a part span shrouded titanium fan blade 
was subjected to vibration testing. The test objectives were to measure the ma- 
terial damping in the blade for its first three resonant modes and also to assess 
the effects of elevated temperature on the damping characteristics of the blade. 
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3 . 3 . 3 . 1 Test Approach 

The analysis of Section 3.1 requires the determination of the principal 
stresses throughout the component. However, if the responses are modal, then 
the stress at a reference point on the blade can be used to scale the modal 
stress distribution. Therefore the test method chosen was to harmonically ex- 
cite the fan blade which was welded at its root to a massive block and measure 
the strain response for various levels of forcing using miniature foil strain 
gages mounted at high axial stress locations near the root and the A.S.M.T. 
position. The loss factors for the first flap mode were estimated from freq- 
uency sweep response plots using the half power method. For the two higher modes 
the loss factors were estimated from strain transient curves resulting from abrupt 
cut-off of the harmonic excitation (i.e., in an essentially single mode response). 
The basic damping data were obtained with the blade vibrating in evacuated con- 
ditions. However, for the tests in a thermal environment, evacuated conditions 
were considered unnecessary since the measurements of any changes in damping 
characteristics due to temperature would provide the desired information. The 
tests were therefore performed in air with temperatures on the blade of the order 
of 150°C. 

The test specimen was an advanced part span shrouded titanium (Ti) fan blade 
with dovetailed root configuration. The same blade was used throughout the program 
to obtain measurements of damping at root as well as at shroud interfaces. 

For this test series, the blade was welded into a massive titanium block at 
its root by machining the dovetail off and replacing it with weldment (see Figs. 

16 and 17) . The root was oriented in the block such that the applied excitation 
vector was normal to the tip chord. The system was mounted in the f ixture/vacuum 
chamber vibration set up described in Appendix B. 

The test arrangement for the thermal testing is shown in Fig. 18. It was 
basically the same as that used for the "cold” blade except that the airtight 
alumimum lid was replaced by a fiberglass cover and two 1 kW blowdrier heater 
elements were installed to provide radiation heating. The heaters were freely 
suspended from a separate bracket arrangement such that they were approximately 
215 cm from the blade profile. Three thermocouples were used, one inside the 
dovetail of the blade, one near the shroud, and one close to the tip. A VARIAC 
voltage controller was used in all the tests and a prolonged soak (three hours) 
assured a stabilizing of the temperature profile from root to tip. The two TCs 
on the blade touched the blade surface before the tests started so that the temp- 
eratures recorded were actual metal temperatures, but were moved away during vi- 
bration. 

3 . 3 . 3 . 2 Instrumentation 

The blade was instrumented with two high temperature strain gages (BLH 
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Electronics, Type FSM 12-:33-.S6) . The lead wires were orignially attached using 
M BOND GlOO high temperature cement which proved to be too brittle and was re- 
placed with RTV 106. The choice of gage location was directed by the need for 
high axial strains in all three modes and a NASTRAN dynamic analysis was used to 
determine the optimum locations. These locations, as shown in Fig. 16 and 17, 
were near the root at the leading edge and above the shroud near the maximum 
thickness position. 

The instrumentation system used in the tests was the basic system described 
in Appendix B. An amplifier was required in the input circuit to enable the low 
(< O.lg) acceleration levels to be applied. The input control accelerometer was 
mounted outside the fixture at the root location (see Fig. 18). For the tests, 
all output voltage RMS readings were taken from the digital voltmeter and these 
levels were noted on the data charts. 

3. 3. 3. 3 Tes ting and R esult s 

Testing at ambient temperature consisted of exciting the blade at each of its 
first three natural modes of vibration under evacuated conditions down to 10 torr 
and recording the strain response. Decay tests in the first mode (first flap) 
were not feasible owing to difficulties experienced in achieving a total arrest 
of input motion upon abrupt cut-off of the excitation. This difficulty was caused 
by the inertial forces of the vibrating blade being sufficiently large to induce 
vibratory motion in the shaker armature after the input was cut-off. Thus a mod- 
ified procedure using a frequency sweep over the resonance and determination of 
the bandwidth at the "half power point" at various degrees of pressurization and 
evacuation was utilized. The pressure was gradually reduced until the damping 
limit imposed by the available minimum sweep rate (0.0136 octave/minute) was rea- 
ched. This damping value was calculated to be q = 0.0022. The loss factor for 
material damping, being much smaller, was obtained by extrapolation of the data 
points to zero pressure. The values obtained for the first mode loss factor and 
frequency are shown plotted against the pressure ratio in Figs. 19 and 20 respec- 
tively. The loss factor value obtained by extrapolation can be seen to be between 
0.003 and 0.00085 for a range of "root" stresses of 56 to 110 MPa. These results, 
together with the data for the two higher modes are tabulated in Table III. It 
should be noted that the loss factors for modes 2 and 3 were obtained from the 
average slope between points 1 and 4 dB down from the maximum value on the decay 
curves. Figure 21 shows a comparison of the test results (loss factor) with the 
theoretical values for the first three modes. 

In the thermal environment, the loss factors in the first mode were high 
enough, because of air damping, to be derived directly from the slow frequency 
sweep response data. Results for the first mode are shown in Fig, 22, Two 
input levels, namely 0.075 g and 0.15 g, were applied and loss factore were 
determined using the "half power method," As can be seen, the temperatures 
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dropped during the test as a result of air being pumped around inside the 
chamber by the blades and being cooled by the surfaces of the chamber. "Hot" 
and "cold" peak decay curves for the second and third modes are shown superim- 
posed for comparison in Fig. 23 and 24. 

3.3.4 Composite Material Damping 

In addition to the several mechanisms of damping discussed so far, two 
additional mechanisms, viz, (1) composite materials and (2) coatings typically 
used in turbine blading, were examined to determine the potential, if any, of 
damping available. The effort related to composite materials is discussed below. 

Potentials of damping in composite material were evaluated by tests conducted 
on two materials: ' Boron/Aluminum (50% by volume of Boron filaments in a 6061 
Aluminum matrix); and Boron/ Aluminum/Titanium (0.076 mm thick titanium cladding 
on two sides amounting to approximately 6% by volume titanium) . The specimens 
were flat strips (see Fig. 25) of nominal size of 20.32 x 2.54 x 0,25 cm. Two 
of each specimens were tested as cantilever beams (see Fig, 26) and data were 
obtained in their first two modes. All tests were conducted in a vacuum chamber, 
exciting the beams at their resonances and recording the decay of vibrations 
upon abrupt cutoff of input excitation. The strain signals from a gage mounted 
0.64 cm from the root were filtered before being recorded on a level recorder. 

Before the final sets of data were recorded to obtain an estimate of 
damping, it was found necessary to conduct tests to (1) determine the clamping 
torque above which the fixture participation could be considered to be minimum 
and (2) establish the level of maximum strain below which the material behavior 
is linear. Based on this series of tests, a clamping torque of 60 Nm was found 
to be adequate and about 500 micro strain at the root of the cantilevered spec- 
imen was found to be the desired limit. During these initial tests the material 
characteristics of the specimens appeared to undergo continuous change and this 
necessitated subjecting one B/Al test piece and one B/Al/T test piece to pro- 
longed testing over a wide range of input levels. In doing so, it was suspected 
that these two beams were overstressed and thus were not used to obtain the final 
results. The results presented in Fig. 27 were obtained for bending strains 
within the linear range. 

Using the carefully measured values of resonant frequencies and from a knowl- 
edge of the beam characteristic numbers, the modulus of elasticity for B/Al and 
B/Al/Ti specimens were calculated to be 199 GPa and 166 GPa, respectively. 

Figures 27a and 27b show the measured loss factors as a function of maximum root 
strain. Factors obtained from characteristic beam functions (Ref. 48) were used 
to estimate the root strain from the measured strain. For both the specimens, 
the loss factor in mode 2 was found to be smaller than that in mode 1. The 
variation of damping with strain level is small indicating that the value of 
n is the relationship D = Ja^ is close to 2. Even with limited maximum strain 
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levels, repeatability, particularly in the first mode, was found to be difficult 
as is evident from a comparison of results for mode 1 obtained before and after 
measuring mode 2 . 

3.3.5 Tests on Turbine Blades With Thermal Barrier Coatings 

The objective of this series of tests was to evaluate the damping potential, 
if any, of some typical coatings currently in use on turbine blades. In general, 
these coatings are used for thermal protection and to inhibit corrosion. The 
coatings are applied to the blade airfoil and upper surface of the platform. 

3 . 3 . 5 . 1 Tes t Appr oach 

The approach was to measure the loss factors of blades (before and after 
applying a coating) using response decay techniques. Two turbine blades, sim- 
ilar to that used for the platform damping tests (Section 6 ) were prepared and 
tested. The following treatments were used in preparing the blades for these 
tests . 

Treatment // 1 - Two blades 

Strip (20% HNO 3 ) and vapor blast 
Treatment // 2 - Blade //I 

Diffusion Aluminide coat (pack cementation process) to a thickness of 
76 microns. 

Diffusion heat treat at 1079°C/4 hrs/Argon 


Treatment // 3 - Blade #2 

Sputter coat with NiCoCrAlY to a thickness of 126 microns. 
Diffusion heat at 1079°C/4 hrs/Argon 

Treatment // 4 - Both Blades 


Overcoat (plasma spray) with Y 203 ~Zr 02 to a thickness of 250 microns 
using a 76 microns thick bond coat of NiCr/Al composite. 


Before each test, the blades were instrumented with a strain gage (MM Type 
EA-06-062AK-120) at the root of the shank. 1 

Each blade was tested in the same fixture as was used in the platform damping 
tests (see Fig. 4), but in these tests, additional care was taken to make the 
fixture participation an absolute minimum and torquing procedures were standardized. 
Even with this additional effort, the fixture participation and root effects were 
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still significant. However, initial tests indicated that small changes in blade 
damping could be detected with the setup as long as the torquing procedures were 
followed . 


Strain gage signals were conditioned using the strain gage translator and 
were recorded on the level recorder. Calibration of the level recorder was 
performed by exciting the blade with the shaker and reading the output levels 
on a digital voltmeter. 

Decays were obtained by impacting the blade with a 'soft’ hammer rather 
than by using the decay from resonance technique as tuning was impractical 
with the extremely low damping levels prevailing. 


3 . 3 . 5 . 2 Jfes t Resud. t s 

The ranges of loss factor and response frequency for the blades following 
a treatment are shown below for a stress range of 17 MPa - 55 MPa. 

Blade //I Blade #2 

Treatment # Frequency Loss Factor Frequency Loss Factor 

Hz Hz 

520.17-520.09 . 00030-. 00031 523.25 . 00029-. 00037 

521.16 .00022-. 00029 

521.27 .00024-. 00036 

.00090-. 00128 501.76 .00060-, 00106 


1 

2 

3 

4 


3.4 Numerical Results and Discussion 

In this section, numerical results relevant to the evaluation of material 
damping in titanium are presented. Theoretical results are compared to those 
obtained from tests. 

3,4.1 Damping Constants 

The measured data from tests on the uniform cant;ilever beam are given in 
Table I. The dimensions of the beam are given in Fig. 7a. The following prop- 
erties have ^een assumed for titanium, the beam material: E = 1.27 x 105 MPa; 

p = 4370kg/m, o - 590 MPa. Use of the test data (excitation vs response) in the 
computer program DMPCON yields the following results: 
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Mode 

Frequency 

Theory 

(Hz) 

Test 

n 

3 

kN.m/m /Cycle 

1 

67.1 

68.4 

1.89 

1.30 

2 

420 

427 

2.05 

2.40 

3 

1178 

1196 

2.02 

1.03 


Average Values ~ 2 1.58 

It is noted that n is very close to 2 for titanium. 


An analysis of data from tests on the uniform cantilever beam indicates 
that the value of the damping constant n is very close to 2, for titanium, 
suggesting that material damping is essentially independent of stress. This 
value of n is in general agreement with the value quoted in Ref. 2. However, 
the value of J* = 1.58 kN.m/m /cycle that was obtained in this program, is an 
order of magnitude less than that noted in Ref. 2 (J* = 14). This considerable 
difference could be attributed, among other factors, to differences in test 
frequencies. It is noted that values of damping constants in Ref. 2 were obtained 
from tests conducted at 0.33 Hz, essentially static conditions. It is necessary, 
however, to emphasize that no final conclusions should be drawn from this limited 
study. Only a thorough program solely devoted to a study of material damping 
constants for a set of selected engine blade materials can serve to establish 
the magnitudes of these constants for use in design considerations. 

3.4.2 Twisted Plate 


As indicated in Section 3.1, the procedure for calculating material damping 
calls for the modal stress distribution of the component. The NASTRAN finite 
element model of the twisted plate is shown in Fig. 10. The model consists of 
48 QUAD4 type of elements. Figure 28 shows, for the first three modes, the 
spanwise distribution of modal principal stresses near the midchord of the 
plate. A comparison of the measured and theoretical (NASTRAN) frequencies (Hz) 
is shown below: 


Mode 

Zero 

Theory 

Speed 

Test 

500 

Theory 

rpm 

Test 

1000 

Theory 

rpm 

Test 

1500 

Theory 

rpm 

Test 

1 

97.8 

101 

100 

103 

107 

109 

118 

119 

2 

540 

550 

542 

553 

548 

558 

557 

568 

3 

670 

667 

670 

677 

671 

677 

672 

678 
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A theoretical loss factor of 0.00018 is obtained by using J* = 1,58 and 
n = 2. Figure 29 shows the test results of the forced resonant response of 
the twisted plate. It is possible to estimate the extent of material damping 
by using the forced response results; what is required is an equation relating 
the material damping loss factor, blade stress, and the base excitation level. 

Such an equation has been derived in Ref. 42 for the case of a uniform cantil- 
ever beam. The necessary analysis involves the computation of the work done by 
the shear force at the root and the strain energy stored during one cycle of 
vibration. For the case of the twisted plate these quantities were calculated 
by using modal information (frequency, mode shape and modal stress) obtained from 
a NASTRAN analysis. An analysis similar to the one in Ref. 42 was carried out 
and the results from Fig. 29 were used to compute the material damping loss 
factors as follows: 0.0001 for the first mode and 0.00016 for the second mode. 

The results of calculating the damping constants by using the uniform 
cantilever beam tests indicate that n is very close to 2 for titanium. This 
implies that material damping for any component made of titanium is essentially 
independent of stress; suggesting that the detailed procedure developed in 
Section 3.1 is essentially redundant if the blade material is titanium. However, 
in order to study the sensitivity of the theoretical results to possible variation 
in the value of n, and to checkout the computer program, MATDMP, loss factors 
for the twisted plate were calculated by using the modal stress distributions 
obtained from NASTRAN. The results for a reference root stress of 68.49 MPa 
are summarized below. 


Mode 

Theoretical 

Loss Factor (p) ; 

J* = 1.58 

n = 1.9 

n 2.0 

n = 2.: 

1 

.00024 

.00018 

.00014 

2 

.00024 

.00018 

.00014 

3 

.00023 

.00018 

.00015 


For the twisted plate, centrifugal forces caused an apparent increase in 
material damping in the first two modes (see Fig. 12) possibly due to some 
participation by the fixture. No such increase was observed in the third mode 
(torsion) indicating a more complete restraint of the fixture against torsion. 
Data were taken at vibratory stress levels not exceeding 46 MPa. Attempts at 
testing at higher stress levels were thwarted by tip deflection limitations 
(driver proximity) in the first mode and driver power limitations for the higher 
modes. This small stress range does not enable the nature of dependence of 
damping on stress to be established. 
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3.4.3 Fan Blade 


A NASTRAN finite element analysis was carried out to determine the natural 
frequencies, mode shapes, and modal stress distribution for the first three 
modes of a typical fan blade. A comparison of theoretical and test frequencies 
is shown below; 


Frequency (Hz) 

Mode 1 

Mode 2 

Mode 3 

Theory 

77.3 

242 

497 

Test 

77.9 

250 

482 


The test results for material damping (loss factors) for the fan blade are 
shown in Fig. 21. A theoretical loss factor of 0.00018 is obtained by using 
J* = 1.58 and n = 2 for titanium. 

For the fan blade, test results in the first mode indicate some spread in 
data, possibly due to the method that had to be used to obtain these results; 
i.e., the method of measuring damping by vibrating the blade in a pressurized 
chamber and extrapolating the results to vacuum conditions. Measurements in 
modes 2 (second bending) and 3 (first torsion) were obtained more directly by 
observing the decay of vibrations in vacuum. 

Tests performed on the fan blade in a thermal environment indicate that, 
for titanium, the differences in material damping between "cold" and "hot" 
conditions are insignificant. This appears to confirm the results of Ref. 12. 
The observed drop in resonant frequency in each mode, with increase in temper- 
atures, is attributable to the decrease in modulus of elasticity at higher temp- 
eratures . 

B/AL and B/AL/Ti specimens tested show loss factors between .001 and .002, 
although repeatability of tests was a major problem. Tests performed on turbine 
blades coated with typical corrosion inhibitor coatings show no change in 
damping levels. However, for blades having a 250 microns thick thermal coating, 
an increase in damping of 300 to 400 percent was noted. 
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4. DAMPING DUE TO SLIPPING AT ROOT 


The root structure of a jet engine blade is typically of a dovetail, 
pin or firtree design. The root slides into a corresponding slot on the rotor 
disk, and centrifugal forces provide a tight fit during engine operation. 

Major parameters which control the extent of friction damping at the root 
to disk interface are: centrifugal loading, roughness of surfaces, level 

of external excitation, and root geometry. Studies which address the 
specific problem of root damping in a jet engine blade can be found in 
Refs. 18-22. Presented in this Chapter, are the results of analysis and 
testing of a fan blade having a dovetail root. 


4.1 Forced Response of a Fan Blade with Slip at the Root 


Based on an assessment of different available analytical approaches, 
the approach of Ref. 20 was chosen for detailed study and application. The 
analysis presented in this Section is based on the following assumptions; 
see also remarks in Section 2.1.2: 

1. Friction mechanism is of macroslip-type, and governed by Coulomb’s 
law of dry friction. 

2. Friction damping does not drastically alter blade resonant 
frequency and mode shape, so that attention can be focussed on 
a single chosen mode of vibration. 

3. For the mode under consideration, the modal stiffness and modal 
stress are normalized to unit tip displacement of the blade, and 
are computed for the case when there is no slip at the root. 

For a vibrating blade with slip at the dovetail root, the total 
motion may be considered as a combination of rigid body motion and elastic 
vibratory motion as shown in Fig. 30a. If attention is focused on a single 
mode of vibration, then the governing equation of vibratory motion can be 
written as 

rri| X, + k(x,-X 2 ) + c(x, - X 2 ) " ^ • (4.1) 

Here, mj^, k, and c are modal quantities, P is the external excitation 
amplitude, co is the excitation frequency, x-]_ is the tip amplitude, and X 2 
accounts for blade rigid body motion. For the rigid body motion the 
equation is 
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IqB +Mf sgn (e)+My = 0 , 


(4.2) 


where Iq is the mass moment of inertia of the blade about the root, Mf 
is the frictional moment due to slip, and is the moment due to 
vibratory forces. Letting, 


Iq “ ro2 I- , 

Mf = R/l N cosa , 


Xg " » 

K = (r/l) cos a, 


(4.3) 


where L is the length of the blade, y is the coefficient of friction, N is 
the axial pull on the blade, a is the dovetail angle, and R is the "effective 
radius" to the rubbing surface (Fig. 30b), the equations of motion for the 
blade are 


m|X|+ k(X|-X2) + c(x,-X2) = P cos cut , 
m 2 ^ 2 “ k(x,-X2) -c (x, -X2) + /X N sgn(x2) = 0 , 


( 4 . 4 ) 


where overbars denote dimensional quantities. It is convenient to non- 
dimensionalize the equations of motion by introducing the following 
quantities 


^ __ 2 “ — 

/3 = m2 /m| , >^i = X| /xq , X2= X2/X0 , ojo "" k /m, , t = tojQ , 


OJ 


= to /djQ , F=/iKN/kxQ,p = P/kXo, C = c/ 2 m,coo , 


(4.5) 


where Xq is some characteristic length and the governing equations in 
dimensionless form are. 


X| + (X1-X2) + 2 ^ (x, -Xg ) = P COS cut, 
jSx 2 ~(xi-X 2 )~ 2 ^(X|-X 2 )+ F sgn (X2) = 0 . 


(4.6) 


Assuming a solution of the form 

X, = a COS (cot + - 8 ), 

X2 = b cos (cot + <^), ( 4 . 7 ) 

sgn(x2)~ “ (4/77-) sin (cot +(/>), 
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leads to 


770,+ {l-o;^)a2“Q2 " cos (f> , 

(a;^- Da, + 7702 - 17Q2 = “Sin(^, 

77a, + a2-(l~/Qaj^)Q2 =0, 
a, - 7702 + 77Q2 + G =0, 

where, 

77 = 2 ^cu, Q2 = b/p, G= 4 F/ 7 Tp, 

0| = (o/p) sin o, 02 = (a/p) cos 8. 


Equations (4.8) are set of algebraic equations in the unknowns, a^, 0^2 , <j>, 
and 6. The last two equations of the set can be shown to yield 


a,= - (17/80)^03+ 6 )/{i+ 7 j*), 

a 2 = [( 1 + 77^ -jSuj^) Q2+ 77 g]/(i + 77^) 

so that , 


(4.10) 


0i = a/p = (af + a| and tan S =a^/a2- ( 4 . 11 ) 

Substitution of Eqs. (4.10) in the first two equations of Eq. (4.8) leads to, 

ton(^ = (oj^Oi + G)/(a2+/9Q2)t*^^> (4.12) 


and after some algebra, a quadratic equation, the solution of which yields. 


where 


Q 2 ~ ^2 ^ 


(CoG)^- c.Cc.g^ 



(4.13) 


C, = 77^(1 +yS)^+ (I+ZS-tSoj^)^, 0^= 77/C, , 

C 3 = [{\~(x)^f + 77^]/oj'^C, , C 4 = (l + 77^)/cj'^C,. 
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It is noted that is a measure of the blade response and Q 2 a measure 
of the slip amplitude at the root. From Eq. (4.12), the condition for slip 
to occur (i.e., Q 2 > 0 ) is 


( 036 ^- 04 ) <0; G < + (4.14) 


When there is no slip at the root, the solution is given by 


Q 0 •" 0 > 


Q, =[(1-0,^)^+.,^]'^^, 


(4.15) 


The analysis and solution developed above is applicable to the case 
of direct excitation of the blade. If the excitation is applied at the 
base of the blade, as was done for tests performed in this program, the 
analysis follows along similar lines to the one given above. The 
necessary modifications to the analysis and solution are as follows: 


P = P/Xo , C 4 = I . 


(4.16) 


a, = -(l+/5-ySo;2)G/C,, 

az- [ 17 (i+/S) 6 /C|] + Qg, 


(4.17) 


where p is the displacement amplitude of the base excitation. 

The influence of friction at the root can be expressed as a loss 
factor, defined as 


' D 
■ 2 7T U 


(4.18) 


where D is the energy dissipated per cycle of vibration and U is the peak 
strain energy in the blade. The two energies are given by 


D = 4^k N X 2 , 



(4.19) 
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so that, in terms of dimensionless quantities 


l7,»6Q2/0f, (4.20) 

and total loss factor for the blade is, 

’^b^Vo'^Vf (4.21) 

where rio accounts for damping from sources other than friction, such as, 
material hysteresis and air damping. 

The major interest is in obtaining the maximum blade response at 
resonance. One approach is to calculate the response for different values 
of the nondimensional frequency, o) and pick out the maximum amplitude. 
However, for the damping levels expected in the problem being studied, 
maximum response occurs when the excitation frequency is near the natural 
frequency (i.e., oa Ri 1) so that a good approximation to the maximum response 
is obtained by simply getting co = 1 in the solution. 

The blade tip deflection can be related to the blade stress by 
specifying a baseline or reference condition where the blade vibration 
characteristics are fully known. If the baseline tip deflection of the 
blade is chosen as the characteristic length x^ used for nondimensionaliza- 
tion in Eq. (4.5), then it can be shown that 


cr - Q Cq , 


(4.22) 


where is the blade stress in the baseline condition and a is the blade 
stress at any other condition. In the present study, the condition of no 
slip at the root is chosen as the baseline condition so that the baseline 
balde stress is the same as the modal blade stress, 

A computer program, ROTDMP was written to compute blade loss factor 
and blade stress (near the root) as a function of blade axial load. The 
program calculates the resonant response by using Eqs. (4.10) through (4.13), 
the loss factor by Eqs. (4.18) through (4.21), and the blade stress by 
Eq. (4.22). The input to the program includes: frequency in Hz; blade 

axial load, N; coefficient of friction y; dovetail correction factor k; 
mass ratio 3; modal stiffness k; baseline (modal) stress Oq and non-friction 
damping loss factor tiq. 
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4.2 Experimental Investigation of Damping at Blade Root 

The objectives of this series of tests were to estimate the potentials 
of damping available in a root structure of dovetail design and to compare 
the results with the trends established by analysis. 

4.2.1 Experimental Approach 

Section 4.1 describes the analysis of the forced response of a single 
fan blade with slipping in a dovetail configured root. Since the major 
effect of the centrifugal loading of a typical blade is to increase the 
reactive forces on the mating surfaces of the dovetail, the approach taken 
in this investigation was to test a single dovetailed blade and simulate, 
as closely as possible, the root loading by means of a special fixture. 
Sinusoidal excitation was applied at the blade root and the blade response 
was monitored via strain gages mounted at high stress points. Since the 
major interest was in obtaining the maximum blade response at resonance, 
the tests were run at steady state at a resonant frequency and transient 
decays were obtained by abruptly cutting off the excitation. The level of 
damping was estimated from the slope of the resulting decay curve. 

The test piece used was another specimen of the titanium fan blade used 
in the material damping tests. The blade was mounted in the fixture, 
described in the next section, using two titanium alloy (8-1-1) clamps with 
mating surfaces having geometry and dimensions exactly the same as those 
of the blade root slot in an actual disc. Blade strain response was measured 
by two strain gages, one mounted near the leading edge at the root and the 
other at the above-shroud-maximum-thickness position. These were the same 
locations used in the material and shroud damping investigations. 

The basic fixture, slip plate and shaker arrangement, as described 
in APPENDIX B, was used. A new steel root block was designed to provide 
a dovetail slot for the blade root and a means of applying an axial load 
to the blade root. The arrangement is shown in Fig. 31. 

The rubbing surfaces were provided by two titanium clamping blocks 
bolted to the steel fixture block. The width dimension of the dovetail 
cutout is maintained by a loading block separating the two clamping blocks. 
The loading block is a sliding fit between them. Strain gaged loading screws 
acting on the block provided a measure of the load along the axis of the 
blade. The load is transferred to the blade root via a small diameter 
roller bearing located on the centerline of the loading block. The assembly 
was then bolted into the vacuum chamber (see Fig. 32), which in turn was 
mounted on the slip plate driven by the electrodynamic shaker. 
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The load bolts were designed to provide adequate axial stiffness at 
the root, yet be sensitive enough to be able to apply accurate loads as 
monitored by means of strain gage output. The steel bolt "working" section 
was 2.54 cm long and 0.32 cm in diameter. Two strain-gages were used on 
each bolt and arranged to minimize bending effects. As a result of 
calibration tests, .695 microstrain per Newton was determined. 

4.2.2 Instrumentation 


The same type of strain gages as previously used in the material 
damping tests were used on this blade. The instrumentation system used 
was the same as that discussed in APPENDIX B except that the X-Y recorder 
was not required. All strain gage responses were recorded as unfiltered 
and filtered signals. The Krohnhite filter was set at 100 Hz in the low 
pass mode to obtain the filtered results. 

The two strain gaged loading bolts were connected through a 3 pole 
switch to a portable balance bridge and strain indicator unit as shown 
in Fig. 33. 

4.2.3 Test Results 

All tests were performed with the pressure in the chamber reduced to 
15 torr. 

After several preliminary tests and developing the technique of 
adjusting the load bolts, the first mode responses were obtained with root 
loading varied from 576 to 2879 Newtons and with excitation levels of 
0.2, 0.5 and 1.0 g applied. The results of the final test series are 
tabulated in Table IV. Typical decay curves are shown in Fig. 34 for the 
minimum root load applied and in Fig. 35 for the maximum root load. The 
loss factors were calculated from the average slope between points 1 and 
4 dB down from maximum on the decay curves. Figures 36, 37 and 38 show 
how the loss factor, root stress and response frequency vary with axial load 
at root and input excitation level. 


4.3 Discussion of Results 

In this section, a comparison of the numerical results from analysis 
and testing of a fan blade with a dovetail root is presented. Values of 
the parameters used in computation are as follows; 
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Blade Length, L 

Mass moment of inertia about root, Iq 

Dovetail angle, a , 

Dovetail effective radius, R . . . , 
Dovetail correction factor, k . . . 

Coefficient of friction, y 

Non-friction damping loss factor, po 

Frequency, (Test) 

Modal mass , m^ 

Modal stiffness, k 

Modal stress near root, 

Mass ratio, 3 


0.28 m (11 in) 

(0.066 slugs-in^/12) 
45 deg. 

8,26 mm (0.325 in) 

0.0209 

0.3, 0.35 

0.0005 

74 Hz 

(2.734E-3 slugs/12) 
11.21 kN/m (64 Ib/in) 
298 MPa (4.35E+4 psi) 
2 


The modal information is for the first mode, and normalized to unit tip 
deflection of the blade. 

The theoretical and test results of blade loss factor as a function 
of blade axial load are shown in Fig. 36 from which the following observations 
can be made. 


(1) Loss factors decrease rapidly with increase in axial load. 

(2) For low axial loads, higher levels of excitation lead to higher 
levels of friction damping. With increased axial loads, this 
dependence on excitation is reduced. 

(3) Loss factors begin to level off at moderate axial loads, asymptoti- 
cally approaching the non-friction damping loss factors; i.e., 

the blade tends to reach the condition of essentially no slip at 
the root. 

(4) At high axial loads, relative macroslip at the root interface is 
unlikely to be present and friction damping, if any, may be 
governed by microslip or local slip. 

Figure 37 shows the theoretical and test results of blade stress 
near the root as a function of blade axial load. It is seen from the figure 
that, for a given axial load, the theoretical blade response (hence blade 
stress) is independent of the excitation level. This is in agreement with 
the results published in Ref. 20. However, test results from the present 
program indicate that, for a given axial load, blade stress is dependent 
on the excitation level. This discrepancy between theoretical and test 
results is possibly due to the idealized nature of the analytical model. 

As the analytical model considers only macroslip motion, an increase in 
excitation level leads to an increase in slip amplitude at the interface. 

The resulting increase in friction damping tends to maintain a constant blade 
stress. The dependence of stress on excitation level, as indicated by test 
results, implies that during tests the slip motion at the interface may 
have been predominantly of the microslip type. 
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5. DAMPING DUE TO RUBBING AT SHROUD INTERFACES 


In the case of part-span shrouded blades, untwisting of the blades under 
centrifugal pull brings neighboring blades into contact at the shrouds. Deter- 
mination of the precise nature of conditions at such an interface continues 
to be an unresolved problem defying both analysis and measurement. Major 
parameters which control the nature and extent of friction damping at the 
shroud to shroud interface are: centrifugal loading, roughness of surfaces, 

level of external excitation, shroud location on the blade, blade geometry, 
and shroud geometry. Studies which address the specific problem of shroud 
damping can be found in Refs. 25-29. Presented in this chapter are the re- 
sults of analysis and testing of a part-span shrouded fan blade with rubbing 
at its shroud faces. 


5.1 Forced Response of a Fan Blade 
With Rubbing at the Shrouds 

Uncertainties in the precise nature of vibratory kinematics and kinetics 
at the shroud faces during rubbing, demanded that different analytical 
approaches be explored to determine their suitability. The microslip approach 
was abandoned in favor of the macroslip approach, after a NASTRAN finite ele- 
ment analysis confirmed the need for an iterative-type solution to obtain 
accurate quantitative results; see remarks in Section 2.1.2, and also Ref. 25. 
Apart from some questions as to its utility in a practical application such an 
iterative-type solution was beyond the scope of the present program in terms 
of time and cost schedule. In the macroslip approach, initial efforts con- 
sisted of modeling the rather high inplane stiffness of the shroud as a spring 
in series with the friction forces at the shroud faces. Such a modeling 
procedure leads to the following result: for a given excitation level, there 

is a rapid decrease in blade stress with increase in shroud normal load. How- 
ever, test results indicated precisely the opposite trend, i.e., an increase 
of blade stress with increase in shroud normal load, for a given excitation; 
thus revealing the inadequacy of the initial analytical model. The difficulty 
is one of proper modeling of the physics of a vibrating blade having rubbing 
at the shrouds. Details of the model that was finally adopted and the assoc- 
iated analysis are given below. 

Schematics of a part span shrouded fan blade, and shroud rubbing surfaces 
are shown in Fig. 39a and Fig. 39b. The following assumptions are made; see 
also remarks in Section 2.1.2, 

1. Friction damping is of macroslip-type, and governed by Coulomb’s 
law of dry friction. 
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2. Friction damping does not drastically alter blade resonant frequency 
and mode shape, so that attention can be focused on a single chosen 
mode of vibration. 

3. For the mode under consideration, the modal stiffness and modal 
stress are normalized to unit tip displacement of the blade, and 

are computed for the case when there is no friction at the interface. 

4. Shrouds are infinitely rigid and shroud load is uniformly distribu- 
ted over the rubbing surfaces. 

5. Slipping motion is resisted by the static stiffness of the blade 
below the shrouds. 

Under the assumptions made, the equations of motion can be written as 
(see Fig. 40) 


mX| + k(x|-X2) + cX| = P cos cu t , 

(5.1) 

^kX2+ 2/iN sgn(x2)”k(x,-X2)=0, 


where overbars denote dimensional quantities, m, k, and c are modal quantities, 
is the blade tip deflection, X 2 denotes the slip at the shroud, P is the 
excitation amplitude, oj is the excitation frequency, y is the coefficient of 
friction, and N is the normal load. The quantity 3k is the static stiffness 
of the blade below the shrouds and resists the slipping motion. Thus, the 
parameter 3 accounts for changes in the stiffness of the system during slip- 
ping motion. It is convenient to nondimens ionalize the equations of motion 
by introducing the following quantities: 


x, = x,/Xo , X2= X2/XQ , (iiQ=k/m, , oj = oj/ojq , t = tojQ, 
£=c/ 2 m,WQ, F = 2 /iN/‘Rxo, P = P/kXo , 


(5.2) 


where x^ is some characteristic length. The equations in dimensionless form 
are 


X| + (x,- X2) + 2^ X,= p COS o;t, 

/5x 2+ F sgn(x2) “ (X 1 -X 2 ) =0 . 


(5.3) 
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Assuming a solution of the form 

X| = Q COS (cut + <^ - 8), 

xg = b cos (ojt + <^), 

sgn (X 2 ) ^ " (4/7t) sin (cut +<^), 

leads to a set of algebraic equations, 

1701+ (i-cu^jag-Qg = cos 4 > , 
(o»^- Dqi + 7702 - 17Q2 ‘ -sin</), 

(l+/S)Q2-a2 =0 

G + a, =0 

where 

T7=2 ^(j, Q2 = b/p, G = 4 F/ 7 TP, 
a, = (o/p)sinS, ag = (a/p) cos 8. 


(5.4) 


(5.5) 


(5.6) 


The last two equations of Eq. (5.5.) yield, 

Q| =~G, 0t2 " ( I"^^)Q 2 , 

so that, 

Q, = a/p = (of + a| and tanS =a,/a 2 . 


(5.7) 


(5.8) 


Substitution of Eq. (5.7) into the first two equations of Eq. (5.5) leads 
to 

tan (jb = [77 (I + / 5 ) Q2 + (I-oj^)g]/ [{/Q + / 3 aj^+ ) Q2 -77G], ( 5 . 9 ) 

and after some algebra, a quadratic equation, the solution of which yields. 


Q2= -CgG +[(C2 g/- CilCjG^-C^)]''^^, 


(5.10) 
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where 

C|= (a)^+/3oi^-/9) + (l + /3)^, C2= Tj/C, , 

C3=[(l-u)^)^+ t;2]/C| , C^M/C, 


It is noted that Q]^ is a measure of the blade response, and Q 2 a measure 
of the slip amplitude at the shroud interface, from Eq. (5.9) the condition 
for slip to occur, i.e., Q 2 > 0, is 


(Cj6^-C4)<Oi G< [d- 0,2)2 + (5.11) 


The analysis and solution developed above is applicable to the case of 
direct excitation of the blade. If the excitation is applied at the base of 
the blade, as was done for tests performed in this program, the analysis 
follows along similar lines to the one given above, the only modifications 
being a redefinition of the nondimens ional excitation amplitude, 

P=P/Xq, and C4 =cjVc, , 

where p is the displacement amplitude of the base excitation. 

The influence of friction at the shrouds can be expressed as a loss 
factor, defined as 

I D 

= (5.12) 

where D is the energy dissipated per cycle of vibration and U is the peak 
strain energy in the blade. The two energies are given by 


D = 8/j.n xa. U = y kxf , 

so that, in terms of dimensionless quantities 


(5.13) 


r}^ = GQ2/of 


(5.14) 
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and total loss factor for the blade is, 


(5.15) 

where accounts for damping from sources other than friction, such as, 
material hysteresis and air damping. 

The major interest is in obtaining the maximum blade response at reson- 
ance. One approach is to calculate the response for different values of the 
nondimens ional frequency, w, and pick out the maximum amplitude. However, for 
the damping levels expected in the problem being studied, maximum response 
occurs when the excitation frequency is near the natural frequency (i.e., 
w ^ 1) so that a good approximation to the maximum response is obtained by 
simply setting w = 1 in the solution. 


The blade tip deflection can be related to the blade stress by specifying 
a baseline or reference condition where the blade vibration characteristics 
are fully known. If the baseline tip deflection of the blade is chosen as 
the characteristic length Xgused for nondimensionation in Eq. (5.2), it can 
be shown that 


O" “ Q > 


(5.16) 


where is the blade stress in the baseline condition and o is the blade 
stress at any other condition. In the present study, the condition of no 
friction at the shroud interfaces is chosen as the baseline condition so 
that the baseline blade stress is the same as the modal blade stress. 


A computer program SHDDMP was written to compute blade loss factor and 
blade stress as a function of base excitation level in g-units. The program 
calculates the resonant response by using Eqs. (5.7) through (5.10), the loss 
factor by Eqs. (5.12) through (5.15), and the blade stress by Eq. (5.16). 

The input to the program includes: frequency in Hz; shroud normal load N; 

coefficient of friction y ; stiffness ratio 3; modal stiffness k; baseline 
(modal) blade stress Qq, and; nonfriction damping loss factor n^. 


5.2 Experimental Investigation of Damping at Shrouds 

The objectives of this series of tests were to examine the nature of the 
response characteristics of a shrouded fan blade and to obtain a measure of 
damping available due to shroud rubbing. 
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5.2.1 Test Approach 


The method described in Section 5.1 enables the calculation of maximum 
response of a single blade subjected to excitation at its base and whose 
shrouds rub on rigid surfaces under a given normal load. The fan blade used 
in these tests was welded at its root to a massive block and provided with a 
means for applying a static normal load on the shroud surfaces. The whole 
assembly was mounted in a vacuum chamber and subjected to harmonic excitation. 
The shroud normal loads and input acceleration levels were varied and the 
stress response at the ASMT (Above—Shroud~Maximum~ Thickness) position was 
monitored. Damping levels were estimated from transient response curves re- 
sulting from abrupt cut off of input signals. 

The fan blade used in these tests was the same titanium blade used in 
material damping tests. The loading platens (Fig. 41) which were used to 
apply normal loads on the shrouds, have a rectangular rubbing surface designed 
to mate with the curved perimeter of the shroud interface. The platens were 
made of titanium alloy 8-1-1 and hardfaced with tungsten carbide to be com- 
patible with the surface treatment on the shroud interfaces. 

The loading disc holding the load platens (Fig. 41) was positioned around 
the blade so that the platens rested on the shroud surfaces. The cable and 
weight system shown in Fig. 42 was used to impose a normal load on the 
shrouds. The cables were positioned 11.4 cm apart on the loading disk. Upon 
reaching the required torque (and hence, the required normal load on the 
shroud interfaces) , the disk was locked in position and the cables were 
removed. Complete air tightness was achieved by capping the access holes on 
the fixture. The above setup assured constant shroud loading through the 
series of tests. 

Actual flexural and torsional stiffness of the blade at the shroud loca- 
tion were measured and were used in the design of the shroud loading mechanism. 
The measurements of these force and moment coefficients were made in a series 
of static loading tests. Moments were applied via rigid extensions glued to 
the shrouds and displacements were measured using dial indicators. The follow- 
ing results were obtained for the several stiffnesses. 




Stiffness 

Measured 

J!X 

Coefficient 

Value 

M.(}) 






K-p 

^6 

114000 N/m 

- ( r 


-e- 

743 Nm/rad 

M, e ^ \ 

1BT 1^ 


5694 Nm/m 


r 

KHe 

412 Nra/rad 

M, M Applied Moments: F Applied Force 


0, (f) Angular Displacement: 6 

Linear displacement 


43 




5.2.2 Instrumentation 


The instrumentation system used was the basic, system shown in Fig. 43, 
and described in Appendix B. All data recording was done on the level 
recorder with steady state strain gage millivolt readings taken directly 
from the digital voltmeter. The response of the blade was obtained from a 
single high temperature strain gage (BLH Electronics Type FSM-12-35-36 , 

3.18 mm) mounted in the above-shroud -maximum- thickness position (see Fig. 17). 
To monitor the shroud motion, two miniature accelerometers (Endevco 2250 and 
Vibrametrics MIOOO) were mounted on a small nonmetallic block which was then 
epoxied to the shroud (see Fig. 44) . The accelerometers monitored the out of 
plane (axial) motion and the in plane motion parallel to the rubbing surfaces. 
Filtering of the accelerometer signals was required at times to assess the 
magnitude of the fundamental component. This was done using the Krolmhite 
variable filter in the low pass mode. A preliminary survey of the blade re- 
sponse was conducted over a wide frequency range, measurements for which were 
made via a miniature accelerometer mounted at the blade tip. 

5.2.3 Test Results 

Several preliminary frequency sweeps were performed to locate the major 
modal responses. The response plot for tip flapwise acceleration of one such 
sweep is shown in Fig. 45. This was performed at a relatively low applied 
torque (5.2 Nm) and so the boundary conditions at the shroud interface repre- 
sented neither a fully locked nor a freely slipping condition. In fact, during 
tests at this low normal load on the shrouds, a rattling type of motion was 
observed at the shroud location. This resulted in a number of peaks occurring 
as a result of the inherently nonlinear nature of the constraint forces at the 
shroud interfaces. However, within the range of frequencies swept, the ob- 
served modes included three predominantly above-shroud modes. Because of the 
low shroud load, these modes occurred at 260 Hz (bending), 475 Hz (second 
bending) and 880 Hz (torsion) . Decay tests were performed from which the 
damping levels in each of these modes were measured. Upon increasing the 
applied torque, it was noted that the number of nonlinear responses diminished. 
Also, the three above-shroud modes manifested themselves (see Fig. 46) at 
278 Hz (bending) , 1110 Hz (torsion) and 1440 Hz (second bending) . It was 
also noted that the response included some participation from the "free" blade 
bending mode at 93 Hz. 

Subsequent testing was focused only on the above-shroud bending and the 
above-shroud torsion modes. These tests were performed with a reduced pres- 
sure of 16 torr inside the chamber. Two methods of obtaining the responses to 
various input levels were used. The first was to employ a slow frequency 
sweep at constant input level for a given shroud load. The frequency range 
was just large enough to at least include the half power points. The sweep 
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rate, was chosen to be low enough not to distort the response curve at the 
existing modal damping levels. A characteristic of these curves was a 
flat peak over a significant bandwidth around the resonance, as shown in 
Fig. 47. The second method was to tune the frequency to obtain the maxi- 
mum response for a given mode. After the maximum strain value had been 
recorded (Digital Voltmeter reading - mv rms) the input was abruptly cut off 
and the decay curve recorded on the peak level recorder. It was noted that 
this latter method produced higher (generally less than 1 dB) maximum blade 
responses than the sweep method. The steady state results were used for com- 
parison with analysis and are tabulated in Table V. Shroud acceleration data 
obtained for certain test points are given in Table V. The loss factor for 
each test point was determined from the decay curves and was obtained from 
the mean slope between the 1 dB and 4 dB points down from the maximum value 
as described in Appendix B. A sampling of the decay curves obtained for mode 
1 is shown in Figs. 48 and 49. Figure 48 shows how the decay pattern changes 
as the input level is increased for a given shroud load and Fig. 49 shows the 
manner in which the decaying signal changes with shroud load and input level. 
Figures 50 and 51 give the ASMT stress as a function of shroud normal load 
for modes 1 and 2. 

After testing, the platens were removed and their rubbing surfaces exam- 
ined. As can be seen in Fig. 52, the majority of wear took place at the outer 
corners of the surfaces. This occurred even though initially the surfaces 
were true and tight against the entire length. It is likely that initially 
some of the new hardfacing may have quickly worn off allowing the load to move 
outward and concentrate on a small area since the loading device was not able 
to automatically adjust to the new contact condition. The distance between 
these contact points was 5.08 cm. However, the fact that the shroud load 
acted on a very narrow region enabled its value to be determined from the 
applied torque. 


5.3 Discussion of Results 

In this section, a comparison of the numerical results from analysis and 
testing of a part-span-shrouded fan blade with rubbing at the shroud faces is 
presented. Values of parameters used in computation are as follows: 


Coefficient of friction, p 0.3, 0.35 

Nonfriction damping loss factor, 0.0005 

Frequency (Test) 280 Hz 

Static stiffness of blade below shrouds 113.8 kN/m 

(650 Ib/in.) 
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Modal stiffness, k 


79.15 kN/m 
(425 Ib/in.) 


Modal stress, ASMT (Above-Shroud-Max. -Thickness) 1003 MPa 

(1.464 X 10^ psi) 

Stiffness ratio, 3 1,44 

It may be recalled that the analysis focuses attention on a single 
chosen mode of vibration of the blade with rubbing at the shroud faces. 

The identification of this mode, for obtaining modal parameters, is made 
difficult by the complex kinematics and kinetics at the shroud rubbing sur- 
faces. Several NASTRAN runs were made to study the effect of different 
boundary conditions on frequencies and mode shapes of the fan blade. The 
runs included the cases of shrouds unrestrained, shrouds fully restrained 
at the shroud faces, and shrouds freely slipping at the shroud faces. Also, 
vibration tests provided insight into the dynamics of shrouds during vibra- 
tion; for example, the mode shape, and the out-of-plane motion of the shrouds 
observed suggested that the second bending mode of the blade with shrouds 
unrestrained is the appropriate mode of vibration to model. Therefore, 
values of modal stiffness and modal stress from this mode were used in the 
analysis . 

The following observations are made upon a close examination of Figs. 45 
through 52: 

(1) The modes of interest correspond to those in which the motion is 
essentially above the shroud (see Figs. 45 and 46). 

(2) The stress amplitude at slip resonance (see flat top in Fig. 47) is 
proportional to shroud load (N) and excitation level (g) . Off resonance, 
however, the stress amplitude is proportional to excitation level only. An 
increase in shroud load at a given g level, increases the frequency and re- 
sponse amplitude at which macroslipping would begin. For a given shroud 
load, an increased g level hastens the slipping action and leads to a wider 
bandwidth of frequencies in which damping can be effective. At a given fre- 
quency and g level, an increase in shroud load leads to an increase in stress 
level . 

(3) Reference to Figs. 48 and 49 suggests that slipping increases with 
g levels as indicated by the progressively steeper slopes at the start of the 
decay. However, the decay curves superpose at lower stress levels indicating 
damping levels which are essentially independent of initial conditions. 
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(4) Reference to Fig. 50 suggests that, for a given shroud normal load, 
the blade stress approaches a constant asymptotic value with increase in exci- 
tation level; i.e., distinct macroslip conditions. The analysis predicts an 
abrupt transition, at low g levels, from a region of no friction damping to 
the macroslip region. The test results indicate a smoother transition at 
about 1 g for 150 N, about 2 g for 300 N, and about 3 g for 400 N; the smooth- 
ness of the transition indicating a region of partial microslip conditions. 

Analysis suggests that blade stresses are sensitive to the coefficient of 
friction; higher the coefficient of friction, higher the stresses at which 
transition to macroslip occurs. The ability of the analytical model to pre- 
dict blade stress appears to improve with reduction in shroud normal load 
and/or increase in excitation level, i.e., conditions representing macroslip 
motion. 

Measurements made in the above-shroud-torsion mode indicate damping levels 
substantially lower than those in the above-shroud-bending mode. The vibra- 
tory motion in the torsion mode causes a slipping motion (rotational) that is 
different from that caused in a bending mode (translational). In the tests 
performed, the rubbing action at the interfaces in the torsion mode may have 
been microscopic as no trend to reach a constant stress level could be observed 
as g levels increased up to 4 g (compare Figs. 50 and 51). 

It would appear that in these tests, dry friction damping was more effec- 
tive in the bending mode than in the torsion mode. This is most likely due 
to the restraint against torsional motion in the present setup being different 
from that offered by neighboring blades in an assembly. 

The wear pattern shown in Fig. 52 suggests that the shroud contact during 
vibration occurs over a very narrow region and even this minimal contact is 
sufficient to cause changes in mode shapes as well as to provide adequate 
damping . 


5.4 Characteristics of Sliding Friction 

Observations made in the shroud damping tests indicate that slip occurred 
between the mating surfaces of the shrouds. Accelerometers located on the 
shroud monitored the shroud motion. 

The analysis of Section 5.1 employs the normal coulomb damping model 
for dry friction damping at the shrouds with an estimated equivalent coeffi- 
cient of friction of 0.3 to 0.35. The purpose of this series of tests was to 
examine the friction forc^-slip relationship due to the surface contact inter- 
action of two test pieces representative of the shrouds of a typical fan blade. 
From this investigation, the value of equival.ent sliding friction coefficient 
was to be derived. 
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5.4.1 Test Approach 


The approach chosen for this investigation was to provide sinusoidal 
relative motion between two test pieces which were held together with a 
constant normal force and to measure the resulting frictional forces directly. 
The range of normal loads, frequencies, and relative velocities to be im- 
posed were to be consistent with those pertaining to the shroud damping in- 
vestigation (see Section 5.2.3). 

The test pieces chosen were a pair of load platen pieces as used in the 
shroud damping investigation. These were rectangular prisms (22.35 x 32.00 
X 6.35 mm) made from titanium (Ti) 8.1.1 alloy with the rubbing face (6.35 
X 22.35 mm) flame sprayed with tungsten carbide (see Fig. 52). One piece 
was gripped with its rubbing surface horizontal in a special clamp mounted 
on a slip plate attached to an electrodynamic shaker such that motion could 
be applied parallel to the long side of the rubbing surface (see Fig, 53). 

The other piece (shown in Fig. 54 as stationary upper test piece) was mounted 
at the end of a pivoted area attached to the shaker support frame (see Figs. 

53 and 54) . The upper test piece was supported longitudinally by two strain 
gaged bolts (previously used in the root damping investigation) and the nor- 
mal load was applied via roller bearings along the top surface of the test 
piece. Lateral stabilization was achieved by locating the tips of the load 
bolts in shallow identations in the fore and aft faces of the test piece. The 
strain gaged bolts were preloaded to 575 N. The sensitivity of the load 
measuring arrangement was 347 microstrain/kN. The input motion of the lower 
test piece was monitored by an accelerometer mounted on the clamp. The normal 
load was applied by putting the required weights on the loading platform 
as shown in Fig. 53. 

Testing consisted of setting the required frequency and input accelera- 
tion and recording the frictional force-slip loop on an oscilloscope for a 
range of normal loads. 

5.4.2 Instrumentation 

The excitation control system used is described in Appendix B along with a 
50 Hz bandwidth tracking filter to condition the control accelerometer signal. 
The strain gages on one of the load bolts were connected to a bridge amplifier 
with internal calibration capability. The filtered input acceleration control 
signal was connected by the horizontal axis of an oscilloscope and the output 
of the bridge amplifier was connected to the vertical axis of the oscilloscope. 
The horizontal axis was calibrated by putting in a known acceleration and the 
vertical axis was calibrated using the bridge amplifier integral calibrator. 
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5.4.3 Testing and Results 


The nominal condition for testing was defined as a displacement of 
0.127 mm DA at a frequency of 280 Hz with normal loads ranging from about 
250 to 700 N. Three loads were used, based on the weights available and the 
dead load of the upper arm; these were 271, 492, and 672 N. The displacement 
was varied from 0.056 to 0.132 mm DA and frequencies of 80, 140, 280 and 
420 Hz were input. 

Initially, the time varying signals from the two accelerometers and 
strain gages were examined on the scope. A composite picture of a typical 
set of traces for one condition is given in Fig. 55. The magnitude of the 
noise in the quiescent system can be seen in the uppermost trace. As can be 
seen, the upper test piece acceleration has a high frequency component. This 
high frequency may be due to the surface rubbing or possibly to impacting of 
loose roller bearings in the upper arm. The component usually accompanied 
the results using the lowest normal load. Using a 250 Hz low pass filter, 
this component was removed and the magnitude of the fundamental can be seen 
compared with that of the input signal. A friction force-slip loop is shown 
superimposed for comparison. The results for all tests performed are given 
in Table VI. 

The equivalent friction coefficient was determined from the measured 
area of the recorded (photograph of oscilloscope picture) loop and a knowledge 
of the input frequency and displacement. The resulting values for friction 
coefficient are plotted against normal load, input frequency, and relative 
maximum velocity in Fig. 56. 
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6. DAMPING DUE TO PLATFORM-DAMPERS 


Platform-type dampers are frequently used in gas turbine blade design 
to dissipate energy through friction generated at the interface between the 
damper and blade platform. The type of damper considered in this program is 
the so-called blade-to-ground damper (Fig. 4). During engine operation, cen- 
trifugal forces bring one end of the damper into contact with the blade plat- 
form, the other end being attached to a relatively static component such as a 
cover plate. Major parameters which control the extent of friction damping at 
the damper to platform interface are: contact load, roughness of surfaces, 

level of external excitation, location of platform on the blade, and stiffness 
of the damper. Studies which address the specific problems of damping at a 
blade platform can be found in Refs. 30-32. Presented in this chapter are the 
results of analysis and testing of a turbine blade subjected to friction damp- 
ing at its platform. 


6,1 Forced Response of a Turbine Blade with a Platform Damper 

Generic models which can be used to analyze blade-to-ground type of plat- 
form dampers can be found in Ref. 16 and 17. The analytical approach used in 
this Section closely follows that of Ref. 30, the generic model being essential- 
ly the same as the one in Ref. 17. The following assumptions are made; see re- 
marks in Section 2.1.2: 

1. Friction damping is of macroslip-type, and governed by Coulomb’s 
law of dry friction; 

2. Friction damping does not drastically alter blade resonant frequency 
and mode shape, so that attention can be focussed on a single 
chosen mode of vibration; 

3. For the mode under consideration, the modal stiffness, and modal 
stress are normalized to unit tip displacement of the blade, and 
are computed for the case when the damper is not in contact with 
the blade platform. 

A schematic sketch of a turbine blade with a platform damper is shovm 
in Fig. 57a, and the associated analytical model in Fig. 57b. The equation 
of motion can be written as 


mx + kx + Cx - P cos o> t - Fjj » 


( 6 . 1 ) 
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where overbars denote dimensional quantities, m, k, and c are modal quantities, 
P is the excitation amplitude, o) is the excitation frequency, x is the blade 
tip displacement, is the modal displacement at the platform, and is the 
force due to the damper. 


Denoting the blade displacement at the platform by q^, the damper spring 
displacement by q^, and introducing 


<#>aV . (6.2) 

it can be shown that 

fx + Ym-Xm . 0 <ST<t 

*^1 'Pd d I - y m , T<3t<TT (6.3) 

y(t) = -y (t+TT). COS T = |-(2ym/5(m), ym’ MN > (6.4) 


where subscript m denotes peak value, N is the damper normal load, and y is 
the coefficient of friction. 

The damper force is a discontinuous function of time, and is shown in 
Fig. 58a as a function of the blade displacement at the platform. Figure 58b 
shows the blade displacem-ent at the platform and the damper spring displacement 
over one cycle of vibration. It is noted that the slip amplitude at the damper 
to platform interface is (2%~ym) • 

It is convenient to nondimensionalize Eqs . (6 . 1) - (6 . 4) by introducing 

X = x/xo , a = Xm/Xo , y = y/xo , b = ym/Xo y cDQ=k/m, oj = a)/a)o, (6.5) 

t = to3o , 6 = k^/k , C = c/2mo5o, p =P/kXo 

where is some characteristic length. The result is 

X + X + 2^x = p COS cut - €y , (6.6) 

^rx + b-a, 0<cut<T 

^ 1 — b , T<ajt<7T 

(6.7) 

y(t) = - y (t + 7 t), cos t = I “(2b/a), b = /xN/c^jjk^XQ ^ 

From the equations developed above, two limiting cases can be identified as 
follows: damper free (no contact with the platform), t= 0, b - 0, and damper 
stuck (no slip at the platform), x == rr, b = a. 
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Due to the nonlinear (discontinuous) nature of the function, y, the 
solution of Eq. (6.6) is obtained by using the method of averaging, a perturba- 
tion method (Ref. 47). The essence of the method is the assumption that, if 
the nonlinearity is small (i.e., g is small), then the amplitude and phase of 
the response may be treated as constants within a cycle of vibration . A solu- 
tion is sought in the form 

X = a(t) COS [t + 0{t)] = 0 COS / , (6.8) 


where the amplitude and phase are assumed to be functions of time. Application 
of the method of averaging transforms Eq. (6.6) to a set of first-order differ- 
ential equations. 


Q = ( I/Ztt) [- 27T^a - e 0 sin^T - tt p sin (^-cut)] , 
Q - (i/27ra)[aT -(a/2) sin 2 t - Trp cos 


(6.9) 


where the case of near resonance, i.e., (jo~ 1, has been assumed. 


Equations (6.9) yield the steady-state response as a set of nonlinear 
algebraic equations. 


-2Ca + R(e,Q,T) = p sin\// , 

- 2Sa + S(e,Q,T) = p cos i//, 


( 6 . 10 ) 


where 

= S = UJ-|, 


( 6 . 11 ) 


R = - (cq/tt) sin^T , S = (€0/27r) ( 2 t- sin 2 t)^ (6.12) 

Though Eqs. (6.10) are nonlinear equations, due to their special structure, it 
is possible to calculate the maximum response by some algebraic manipulations; 
see Ref. 17. Some of the steps are outlined below. 

Squaring and adding Eqs. (6.10) yields a quadratic equation for 6, the 
solution of which is 

8 =(I/2q){s± [p^-(R- 2ta)^]'^^}^ (6.13) 

Since 6 is a frequency parameter, the maximum response occurs when 6 has 
a double root; that is, 

- (R-2£a)^ = 0. (6.14) 
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By using the definitions for t and R given by-Eus. (6.4) and (6.11), respec- 
tively, it can be shown that Eq . (6.14) leads to a quadratic equation for the 
maximum amplitude. The final results are 

= (l/27n7o) {(tt - 4€G) + [(tt- 4eG)^ + 16 e 7r77o G^] | , 

I +(€/47r)(2Tm“ sin Tm)» [l - ( 2 G/QrT^)], 

where is the maximum respose and is the frequency at which the maximum 
occurs, and 

0m= a/p . G = b/p , % = 2^ . 


(6.15) 

(6.16) 


Further, it is seen from Eq. (6.15) that Qm is a function of ho» and G 
which in turn is a function of the normal load on the platform; see Eq. (6.7). 
The optimum value of the damper normal load can be determined by minimii^ing 
Eq . (6.15) with respect to G; that is 

dQpn 


which leads to 

where G* is the optimum value. The associated optimum response and frequency 
can be shown to be 

2G* » cu*=I+(€/4) • (6.18) 

The analysis and solution developed above is applicable to the case of 
direct excitation of the blade. If the excitation is applied at the base of 
the blade, as was done for tests performed in this program, the analysis follows 
along similar lines to the one given above, the only modification being a re- 
definition of the nondimensional excitation amplitude, 

P = p/xo . 

where p is the displacement amplitude of the base excitation. 


The influence of friction at the 
factor, defined as 


I D 

" 2 7T U * 


platform can be expressed as a loss 


(6.19) 


where D is the energy dissipated per cycle of vibration and U is the peak 
strain energy in the blade. The two energies are given by 
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D= |Xm-yml. U = |-kx^, (6.20) 

SO that in terms of dimensionless quantities 

TJf = (4e/ir)G (Om-G)/Q^ , (6.21) 

and total loss factor for the blade is, 

• ( 6 - 22 ) 

where accounts for damping from sources other than friction, such as, 
material hysteresis and air damping. 

The blade tip deflection can be related to the blade stress by specifying 
a baseline or reference condition where the blade characteristics are fully 
known. If the baseline tip deflection of the blade can be chosen as the 
characteristic length used for nondimensionalizat ion in Eq. (6.5), then it 
can be shown that 

(T = a cTq, (6.23) 

where cr^ is the blade stress in the baseline condition and o is the blade stress 
at any other condition. In the present study, the damper free condition (no 
contact between damper and platform) is chosen as the baseline condition. 

It was found, from tests in the damper free condition, that the blade 
shank stress is a nonlinear function of the base excitation level; see Fig. 61. 
In other words, damping due to sources other than platform friction, ri^, is 
a function of the excitation level. It can be shown that, 

T7o “ P°’m/^0 * (6.24) 


where is the modal stress. 

A computer program PLTDMP was written to compute blade loss factor, resonant 
frequency and blade shank stress as a function of the normal load on the platform. 
The program calculates the baseline damping loss factor by using Eq . (6.24), 
resonant response by Eq. (6.15), resonant frequency by Eq. (6.16), loss factor 
by Eqs. (6.19) through (6.22), and blade shank stress by Eq. (6.23). The program 
terminates if the resonant frequency changes from the baseline frequency by a 
specified amount. The input to the program includes: frequency in Hz; base 

excitation level in g-units; coefficient of friction P ; damper stiffness k^j; 
modal displacement at the platform <+)^; modal stiffness k; modal stress Ojj^; base- 
line stress and; allowable frequency shift, expressed as a fraction of the 
baseline frequency. 
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6.2 Platform Damping Tests 


The analysis of Section 6.1 was developed for application to blade-to-ground 
type platform dampers. The forced response of a bladelike structure having 
such a damper is determined in terms of the modal quantities of the structure, 
the normal load applied through the damper, the stiffness of the damper and the 
frictional characteristics of the rubbing interface. In order to verify the 
analysis, a turbine blade with a suitable platform configuration was chosen 
for the test program and a means was devised to provide damping at the platform. 

The test objective was to obtain the frequency response of the blade for 
various levels of input excitation and damper load. The effect of damper stiff- 
ness was also to be determined. 

6.2,1 Test Approach 

In an operating engine, centrifugal forces acting on the damper provide 
the contact load between the damper and the blade platform. In this test ser- 
ies, the contact load was applied and controlled by a lever and weight mechan- 
ism as shown in Fig. 4. The blade was rigidly clamped in a steel fixture and 
excitation was provided at the root by a powerful electrodynamic shaker. This 
method was chosen so that sufficiently high input force was available to over- 
come the large damping forces expected and to obtain realistically high stresses 
in the shank of the blade. The modal response of the blade was measured using 
strain gages located at the blade root on the shank and at the airfoil root on 
the centerline through the shank. 

The test method chosen was to apply a selected damper load and excite 
the blade sinusoidaly at a given input level in the vicinity of the blade’s 
first mode frequency. The frequency was then tuned to produce the maximum 
blade response. 

Since the level of damping produced by the platform damper far exceeds the 
contributions from other sources, the tests were performed in air. The test 
piece was a second stage high pressure turbine blade having a three teeth fir 
tree root configuration and a platform located at approximately 27.7 percent 
blade span as measured from the shank root . The platform had a generous trail- 
ing edge overhang which was convenient for accommodating the damper . The blade 
material was a nickel chromium high temperature alloy. The basic blade geometry 
is given in Fig. 59. 

The damper was a rectangular prism integral with, and cantilevered from 
a T416 stainless steel loading arm. The free end was arranged to rub on the 
underside of the platform, with the load vector normal to the platform and 
damper end face surfaces. Two dampers were made, both with a length of 10.1 mm 
and thickness of 3.2 mm. Their depth measurements were 7.6 mm and 4.2 mm 
giving nominal stiffnesses at the free ends of 10.3 and 48.6 MN/m respectively. 
However, a static test of the smaller damper indicated a stiffness, in the 
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fixture, of 5.5 MN/m. This figure was used in the analysis. The area of 
the rubbing surface was made the same for both dampers, i.e. 2 mm x 2 mm. 

The damper and loading arms were heat treated to HRC 41 to match the hardness 
of the blade. The basic fixturing requirements were that the blade be held 
rigidly at its root with negligible root damping and that it provide a means 
of applying a damper load up to 900N on the blade platform with no free motion 
in the clamping mechanism. The fixture used is shown in Fig. 4 and comprises 
a two piece clamp bolted to an adapter base plate mounted on the shaker. The 
clamp has steps machined in its upper and lower pieces, which contact the 
peaks of the two outer lobes of the blade root fir tree. Four 3/8 inch bolts, 
each torqued to 48 Nm, allowed the blade to be tightly gripped between the two 
clamp blocks. The loading arm of the damper mechanism was pivoted on a pin 
mounted in the lower clamp block. The pivot was designed to be self adjusting 
and no free play was possible. The mechanical advantage of the loading mech- 
anism was 4:1 with the lever loaded horizontally at the top using a cable pul- 
ley and hanging weight system. A counter balance was attached to the base 
plate to prevent lateral loading of the shaker armature. 

6.2.2 Instrumentation 


The blade was instrumented with two strain gages (Micro Measurements type 
EA-06-062AK-120) , one mounted at the shank root and the other at the airfoil 
root as shown in Fig. 59. Both gages were on the center line of the shank. 
Another strain gage of the same type, was mounted diagonally on the side of 
the large damper. 

The instrumentation system used in this test series was the basic system 
described in Appendix B. For the final tests, the output from the blade strain 
gages was read directly from the digital volt meter and tabulated together with 
the input g level and frequency readings. 

6.2.3 Testing and Results 

Prior to vibration testing, the stiffness of the small damper was esti- 
mated with the loading arm mounted on the fixture by statically loading the 
damper tip and measuring the tip displacement using a dial indicator. The 
measured stiffness was 5.5 MN/m. 

During early exploratory vibration testing, the strain gage on the large 
damper was monitored. An example of the traces obtained is shown in Fig. 60. 

It may be noted that at a constant normal load of 89N on the platform, the 
damper strain gage response changes from nearly sinusoidal to one with a much 
more complex wave form as the excitation level is increased from 2.4 g pk. At 
15 MPa shank stress, the signal clearly shows the periodic loss of contact with 
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the platform. At this point acoustical noise was emitted from the blade. At 
4.6 g pk input level, the signal distortion is more pronounced and the noise 
level is increased. 


An initial test was performed to examine the blade alone response char- 
acteristics and to assess the extent of the root damping. The blade response 
to input g levels from 0.06 to 0.9 are tabulated in Table VII and shown plotted 
in Fig. 61. 

Final testing consisted of setting the platform damper load and input accel- 
eration level and exciting the blade in the vicinity of its fundamental (first 
flap) mode. It was found that there was a threshold of both frequency and input 
level below which the blade showed no significant response. Once the threshold 
had been crossed then the blade vibration would continue even with reduction in 
input level back into the threshold region. Thus, this region could not be 
clearly defined but every effort was made to obtain these threshold results in 
a consistent manner. For the highest levels of platform load (>50 N) the reson- 
ant response frequency varied significantly from that of the fundamental mode 
of the lightly loaded blade. This required an increase in the range of frequen- 
cies to include any other resonant modes in the response. The response results 

for the small and large dampers are given in Tables VIII and IX respectively. 

Plots showing the variation of blade stresses with input level for the range of 
platform normal loads are given in Figs. 62 and 63. It should be noted in Figs. 

62 and 63 that the faired curves used to connect the test points are to aid in 

identification of the various load cases and do not represent theoretical results. 
On Fig. 64 the shank stress is shown plotted against the applied platform normal 
load for input values of 0.5g, 0.7g, and l.Og. Test and theoretical values are 
shown for comparison. 

6.3 Discussion of Results 

In this section, a comparison of the numerical results from analysis and 
testing of a turbine blade with a blade-to-ground platform damper are presented. 
Values of the parameters used in computation are as follows: 

0.15, 0.3 
515 Hz 

5472 kN/m (31250 Ib/in) 

245 kN/m (1400 Ib/in) 

0.12 

2293 MPa (3.348E+5 psi) 
see Fig. 61 


Coefficient of friction 

Frequency, (Test) 

Damper stiffness, kd 

Modal stiffness, k 

Modal displacement at platform, 6^ . . . 

Modal shank stress, 

Baseline conditions 
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The modal information is for the first mode and normalized to unit tip 
displacement of the blade. 

Figure 64 shows the theoretical and test results of blade shank stress 
as a function of normal load for three base excitation levels and two values 
of friction coefficient. For the theoretical results shown, the change in 
resonant frequency is less than one percent of the baseline frequency. From 
Fig. 64 the following observations can be made. 

(1) The theoretical results for the case of p= 0.15 are in general 
agreement with the test results. 

(2) Higher excitation levels and/or lower values of friction coefficient 
imply lower rates of stress attenuation, 

(3) For the relatively high stiffness of the damper used in the present 
program, the rate of stress attenuation is essentially linear. 

The influence of damper stiffness on the response characteristics is 
shown in Fig. 65, where the blade shank stress is plotted as a function of 
normal load, for a given excitation level and for four values of damper stiff- 
nesses. For all the results shown in Fig. 65, the change in resonant frequency 
is less than one percent of the baseline frequency. It can be seen that for 
smaller values damper stiffness, there is an optimum normal load where the stress 
is a minimum. For normal loads higher than the optimum, there is a gradual in- 
crease in shank stress. 

Devices such as the platform damper belong to a class of inherently nonlin- 
ear damping mechanisms similar to those of shroud damping. As in the case of 
shroud damping, the vibration characteristics depend on the nature and extent 
of contact between the damper and the blade platform, as well as the excitation 
level imposed on the system. The traces shown in Fig . 60 are a clear indication 

of the changes taking place at the contact region; a nearly sinusoidal response 
changing to one which emitted audible noise. 

Examination of Figs. 62 and 63 suggest that until a certain combination of 
g levels and normal loads are attained, the blade response would be practically 
not measureable. Upon reaching those levels, the response of the blade would 
depend on the magnitude of the normal load. If the latter is "high", the re- 
sponse mode may correspond to a "stuck damper" condition as is evident from 
slightly higher frequencies measured. If the g levels continued to increase, 
then a condition in which the damper breaks contact may be reached. Under these 
conditions the blade "plunges" into a lower mode resulting in the jumps in stress 
as shown in Figs. 62 and 63. 
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A careful examination of Figs. 62 and 63 also reveals that the airfoil 
stresses which begin being lower than that measured on the shank gradually 
become higher for "larger" loads and at higher g levels. This is a clear 
indication of a change in mode shape which points to the need to include more 
than a single mode in the mathematical model. 
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7. GENERAL CONCLUSIONS AND RECOMMENDATIONS 


The series of analytical and experimental programs discussed in this 
report have generated data which forms the basis for understanding the po- 
tential of nonaerodynamic damping in limiting blade vibration levels. This 
investigation begins to address the designer’s need to be able to set speed 
margins of rotors on the basis of predicted damping levels at low order 
resonances. As the exciting forces acting on rotors vibrating at a resonant 
speed are counteracted entirely by damping forces, only an accurate estimate 
of avialable damping can influence the decision in regard to the range of 
acceptable speeds of the machine. Similarly rotor speeds at which blade 
flutter is likely cannot be accurately predicted without a knowledge of the 
contributions from' all sources of damping. 

All test data and analyses in this report are generated for single blades 
and an attempt has been made to outline procedures to extend the results for 
application to an assembly of blades. Undoubtedly future programs must be 
aimed at verifying these procedures and develop a more firm data base that 
will serve as a guideline to the designer. 

Throughout the present program, the tests underscored the importance of 
fixture design. Enormous care had to be exercised to make sure that the fix- 
tures designed and fabricated were rugged enough to withstand the vibration 
levels and to contribute the barest minimum to damping being measured. These 
are tough requirements and cannot always be fully met. Even with the care that 
was exercised in this program, it was clear that the root restraint for the 
twisted plate was found to be less than adequate in the bending modes. 

Measurements of damping of aerofoil type components must always be made 
in vacuum as the contributions from air can lead to order of magnitude incre- 
ments in measured damping. While this requirement sounds simple to implement, 
considerable difficulties in measurement can arise for components with extremely 
low levels of damping. In fact, material damping measurements for the fan 
blade tested in the program required the use of a pressurized chamber. Data 
taken at gradually decreasing pressure levels had to be extrapolated to esti- 
mate damping under complete vacuum conditions. 

Material damping evaluation was made by testing a beam, a twisted plate 
and an advanced fan all made of a titanium alloy. Material damping constants 
(J and n) were determined from test data which indicated a value close to n=2. 
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The latter is in general agreement with earlier reports but the value of 
J found in this program is different from published results by nearly an order 
of magnitude. The initial conclusions are that this difference may be attri- 
buted to the differences in frequencies of vibration between the tests con- 
ducted in this program (50 to nearly 600 Hz) and essentially static condi- 
tions used in earlier reports. These are, of course, preliminary conclusions. 
Only a thorough program solely devoted to a study of damping constants for a 
set of selected engine blade materials can serve to establish the magnitudes 
of these constants for use in design considerations. 

Testing in a limited thermal environment showed that the differences in 
damping between "cold" and "hot" conditions are insignificant. Similar con- 
clusions are drawn for the differences in damping between zero speed and at 
speed. The latter reached nearly 1500 rpm and any changes in damping observed 
is attributed to possible fixture participation. 

The general conclusion (that can be drawn from all the studies on material 
damping conducted in this program) is that damping due to hysteresis in tita- 
nium is very low; the loss factor being of the order of 2 x 10”^ to 5 x 10 
If n = 2, the same values apply for a whole assembly. This suggests that 
damping due entirely to this source should not be accounted for in design 
against flutter and resonant vibration. 

Insofar as damping available at dovetail type root structure is concerned, 
the data shows that rapid reduction in damping loads occur with increasing 
normal load on the root. Thus, in actual practice, for an advanced fan under 
operating conditions the centrifugal pull acting on blades is so high that 
little or no relative motion at the root-disc interaface is likely. Thus, 
contributions from this source of damping at operating speeds is projected 
to be very minimal. Unless design changes incorporating retention techniques 
which permit relative motions at the root location can be examined and intro- 
duced, damping from this source should be viewed to be of little importance 
for design considerations. 

Damping at shrouds and at blade platforms due to rubbing represent by 
far, the most complex mechanism to model analytically in view of the in-, 
herently nonlinear nature of these mechanisms. Testing to measure contri- 
butions from these mechanisms presents its own problems in view of the con- 
tinuously changing boundary conditions at the interfaces as the loading and 
frequency on the device are varied. 

Damping at shroud interfaces appears to be the most likely available 
source of damping for shrouded titanium blades. The level of damping avail- 
able from the rubbing action is difficult to estimate. However, calculations 
made using the analytical model developed in this program indicate that the 
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principal governing parameters are the excitation (g) levels to which the 
blades are subjected and the normal loads acting at the interfaces. Calcu- 
lations show that a loss factor of about .003 can be achieved under a normal 
load of 500 N and a 5 g loading resulting in a vibratory stress of 190 MPa. 
These numbers are valid for a single blade and no extensions to an assembly 
can be made without a knowledge of a modal analysis of a fan assembly vibrat- 
ing in a given nodal diameter pattern. 

Experience in regard to platform damping suggests that the important 
parameters are the damper stiffness and normal load on the damper. Stress 
reductions of nearly 50 percent can be attained with this mechanism. Stress 
attenuation in turbine blades can occur either due to damping and/or due to 
a change of mode shape. The latter arises due to restraint provided by the 
platform damper. 

Under certain combinations of excitation levels and normal load, 
changes in blade mode shapes were observed in this program. Thus, future 
programs in this area must be aimed at including several modes in the mathe- 
matical model. 

Tests conducted on composite material specimens undoubtedly show the 
damping capacities inherent in composites but the principal difficulty in 
testing composites lies in obtaining repeatability and consistency. Only a 
large number of tests can serve to establish reliable information taking 
into consideration scatter levels. 

The program on coating effects was limited to the study of increase in 
damping, if any, that can be obtained by the use of thermal barrier coatings. 
Test data indicated that the differences in damping before and after the 
coating was applied to a turbine blade could be significant but the overall 
damping levels were still very low. 

Although all the tests were conducted on single blades and appropriate 
mathematical models were developed to represent the several mechanisms of 
damping, procedures to apply this data to predict damping levels in an assem- 
bly of blades have been developed and discussed. 
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APPENDIX A 


COMPUTER PROGRAMS 


Numerical results from the analyses developed during the present study 
were obtained by using several simple computer programs written in FORTRAN. 
Following are brief descriptions of the programs, source listings, examples of 
input and examples of output. 


Material Damping 

In material damping studies, two computer programs were used, one for 
calculating damping constants and the other for estimating the extent of 
material damping in a vibrating blade. 

DMPCON 

Vc 

Objective: To determine material damping constants J and n 

Description: The program generates a linear least-squares fit to data obtained 

from tests on a uniform cantilever beam having a rectangular cross 
section. The test data is a table of base excitation g-levels 
vs strain. The integrations needed for evaluating some of the 
constants in Eq. 3.23 are performed in a separate subroutine 
SIMPS by using Simpson’s rule. 

List of Major Variables: 


FORTRAN 

Symbol 

Report 

Symbol 

Units 

Description 

L 

L 

in. 

Length of beam 

H 

h 

in . 

Thickness of beam 

E 

E 

psi ^ 

Modulus of elasticity 

RHO 

P 

slugs/12/in. 

Mass density 

SIGF 

a 

psi 

Fatigue strength 

SGL 


in. 

Strain gage location, distance 




from root 

MODE 



Mode number 

NPT 



Number of test points 

G(I) 



Base excitaion g-level 

EPS(I) 



Strain 

FREQ 


Hz 

Frequency 

N 

n 


Damping constant 

JS 

it 

j 

in. -lb /in. /cycle 

Damping constant 
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Input ; 


L, H, E, RHO, SIGF 

(5E10.4) 

— — — — — — — • ,.i 

SGL, MODE, NPT 

(E10.4, 2110) 

'■ ' — ■ ■ — 

— — . 

G(l), EPS(l), G(2), EPS(2), G(3) , EPS(3) 

(6E10.4) 


M ATDMP 

Objective: To estimate material damping loss factor for a vibrating blade, 

for a given reference blade stress. 

Description: The program uses modal stress distribution of the blade, obtained 

from a finite element program such as NASTRAN. For each element, 
the energy dissipated and the peak strain energy are calculated 
by using Eq . 3.4. Summation of the energies over all the elements 
and using Eq. 3.5 yields the loss factor for the blade. 


List of Major Variables: 


FORTRAN 

Symbol 

Report 

Symbol 

Units 

Description 

E 

E 

psi 

Modulus of elasticity 

SIGF 

^e 

psi 

Fatigue strength 

N 

n 


Damping constant 

JS 


in . -lb / in . ^/cycle 

Damping constant 

SREF 


psi 

Reference blade stress 

NREF 



Element number where SREF is 




defined 

NELE 



Total number of elements 

ID 



Element number 

H 

^i 

in. 

Element thickness 

A 

Ax 

• „ 2 
in . 

Element surface area 

SIGU 

^U 

psi 

Maximum principal stress: upper 

SIGL 

^1 

psi 

Maximum principal stress: lower 

ETA 

n 


Loss factor 


Input : 


E, SIGF, N, JS, SREF, NREF, 

NELE 

(5E10.4, 215) 


ID, H(ID), A(ID), SIGU(ID), 

SIGL(ID) 

(I5,4E15.7) 
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Friction Damping 


In friction damping studies, three programs were used, one for each of 
the following: slipping at the root of a tan blade, rubbing at the shroud 
faces of a fan blade, and rubbing between an external damper and the platform 
of a turbine blade. A list of the major variables used in the three programs, 
ROTDMP, SHDDMP and PLTDMP, is shown below. 

List of Major Variables: ROTDMP, SHDDMP, PLTDMP 


PORTRAIT 

Symbol 

Report 

Symbol 

Units 

Description 

Program 

ROTDMP 

SHDDMP 

PLTDMP 



Hz 

Frequency 

X 

X 

X 




Base excition 







g-level 

X 

X 

X 

MU 

y 


Coefficient of dry 







friction 

X 

X 

X 

KAP 

K 


Factor of dovetail 







geometry 

X 



BETA 

6 


Mass Ratio/ (Stiffness 







Ratio) 

X 

(x) 


K 

k 

Ibs/in . 

Modal stiffness 

X 

X 

X 

KD 

^d 

lbs /in. 

Damper stiffness 



X 

PHID 

^d 


Modal displacement 







at platform 



X 

N 


lbs 

Normal load 


X 

X 

N1,N2 


lbs 

Range of normal loads 

X 


X 

G1,G2 



Pvange of base exci- 

■ 






tation g-levels 


X 


SIGO 


psi 

Baseline blade stress 

X 

X 

X 

SIGM 


psi 

Modal blade stress 



X 

ETO 

% 


Baseline loss factor 

X 

X 

X 

ETA 

% 


Total blase loss fac- 







tor 

X 

X 

X 

EPS 

E 


Parameter 



X 

NPT 



Number of points cal- 







culated 

X 

X 

X 

FTOL 



Tolerance on base- 







line frequency: 



X 




Expressed as a 







fraction 




KODE 



Output Control 







= Print a plot) 







<1 = print only ? 

X 

X 

>: 




(2 = plot only ) 
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ROTDMP 

Objective: To calculate the total blade loss factor and blade stress as 

a function of blade axial load, for a fan blade with slip at 
its dovetail root. 

Description: The program calculates the resonant response of the blade by 

using Eqs. 4.10 through 4.13, the loss factor by Eqs. 4.18 
through 4.21, and the blade stress by Eq. 4.22. 

Input: 



SHDDMP 

Objective: To calculate the total blade loss factor and blade stress as 

functions of base excitation level, for a shrouded fan blade 
with rubbing at the shroud faces. 

Description: The program calculates the resonant response of the blade by 
using Eqs. 5.7 through 5.10, the loss factor by Eqs. 5,12 
through 5.15, and the blade stress by Eq. 5.16. 

Input : 



PLTDMP 

Objective; To calculate the blade loss factor, resonant frequency, and blade 
stress as a function of normal load, for a turbine blade with 
rubbing between the blade platform and an external damper. 

Description: The program calculates baseline damping loss factor by using Eq . 

6.25, the resonant response of the blade Eq. 6.15, loss factor by 
Eqs. 6.19 through 6.23, and blade stress by Eq. 6,24. The program 
terminates if the resonant frequency changes from the baseline fre- 
quency by more than a specified amount. Also the optimum response 
is calculated by using Eq. 6.18. 
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Input : 
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RSS^OAMP ( 1 ) .OMPCON 


LO 

5C 
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— 







— ■ 


q: 







— 











A 







M 


rvj 









•> 







X 


(NJ 







1/1 









0 


<I 

■v» 
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0 


*• 

J- 

>- H- 




0 

■f 




O. QL 




►M 

— 


o 

• 

z z 




1— 

fsj 


ir> 

:r 

» » 




<t 

— 



.-♦ 

UJ UJ 




I— 

1/1 


y/1 

o* 

Q-'Q-' 




H-l 

0 -J 


CL 

D 

O K- 0 >— 




u 

U'v 


UJ 

o 

2:0.x: CL 

m 



X 

— 



• 

»-Z 

• 



u. 

■a fsi 




_J •■-J •“ 

0 




X « 
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13 .MO*-' 

•ft 



0 

o» 


LO 

•> 

1/1 II 1/>II 

-a- 



z 

+ <-* 



Ift 

•>♦-« »♦-* 





— 0 


UD 

in 

u. u. 





NUl. 

1 — 1 — 

•> 

00 

0 — 0 •> 


cr 
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— * 

UJUJ 
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(-• 
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APPENDIX B 


TEST ASSEMBLIES, INSTRUMENTATION AND DATA REDUCTION 


B.l Test Assemblies 

The test program discussed in this report required the design and 
fabrication of a number of test components and fixtures. Table B-1 presents 
a complete listing of the major components of the test assembly used in each 
test series. Design office assembly drawings and parts lists are referenced 
by number for each assembly. Fabrication engineering orders, other than 
those associated with a particular parts list, are also referenced. 

Whenever possible, fixtures were used for more than one test series 
with minor modification. The basic test arrangement used for most of the 
testing comprised a massive root block, onto which the test piece was either 
clamped or welded, bolted into an aluminum plate vacuum chamber and mounted 
on a slip plate so that vibration could be applied horizontally by a Ling 
Dynamic Systems electrodynamic shaker. The slip plate moved on a film of 
oil spread between the plate and the flat horizontal surface of a granite 
block. For the material damping tests on a fan blade, the root block was 
made of titanium alloy 8-1-1 with the blade welded into it (see Figs. 1, 

16-18). A fixture checkout was performed with this test assembly and with the 
blade restrained at its tip. Frequency response sweeps were performed from 60 
to 1000 Hz at input acceleration levels of up to 3 g. These sweeps showed a 
large axial resonance of the fixture at 550 Hz. However, by locating the in- 
put control accelerometer on the root block as shown in Fig. 18, this mode 
could be controlled to give a flat response at both the input location and 
directly adjacent to the blade root from 60 to 1000 Hz. This frequency range 
was adequate for the free blade vibration tests. 

For the shroud damping investigation, a shroud loading disk arrangement, 
(see Fig. 41), was placed around the blade and bolted in position to lugs on 
the fixture side and base plates. Since the torsion mode of the blade with 
shrouds restrained was found to be at about 1100 Hz, additional fixture check- 
out runs were performed up to 2000 Hz. Again, with the control accelerometer 
on the outside of the root block, a flat input level could be maintained at the 
blade root up to 1600 Hz. 

For the root damping tests, a steel root block was inserted in place of 
the titanium block and provision was made to clamp the fan blade at its dove- 
tail with a variable axial load (see Figs. 31 and 32). Response frequencies 
of less than 500 Hz were to be investigated and hence no further fixture checks 
were required. 
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The twisted plate, uniform cantilevered beam and composite beam specimens 
were attached to the steel block individually by special clamps. Because of 
the cantilevered design of these clamps, the input control accelerometer was 
mounted directly on each clamp to ensure correctness of input level. 

The fixturing for the remaining test series is discussed in the body of 
the report in the sections referenced in Table B-1. 

Supporting mechanical equipment to provide the environmental conditions 
and applied loads are included in the list of test equipment given in Table B-2. 


B.2 Instrumentation 

The electronic equipment used for excitation control and for response 
signal conditioning and recording during the test program are included in 
the list of test equipment given in Table B-2. 

The basic control system for all shaker excited tests is shown in 
Fig. 43. The servo level and frequency control was provided by a B & K Auto- 
matic Vibration Exciter Controller which was connected into the shaker and 
matched power amplifier circuit. Test piece input acceleration levels were 
sensed by a VSS Type 101 accelerometer. The accelerometer signal was condi- 
tioned by a charge amplifier located in the LDS Vibration Meter. The output, 
normalized to 10 mv/g, was passed through a 50 Hz bandwidth tracking filter 
provided by the SD Dynamic Analyzer, This ensured that the fundamental com- 
ponent of the input signal had the desired level. For inputting very low g 
levels (< 0.1), a Krohnhite amplifier was interposed between the tracking 
filter and the controller to enable the input signal to be amplified to a 
level acceptable to the controller. The filtered control signal was also amp- 
lified, as required, and connected into the horizontal axis of the oscillo- 
scope in order to obtain Lissajous figures. Logarithmic frequency sweeping 
was provided by the controller with rates down to 0.0136 octave/minute. A 
Fluke counter provided an accurate readout of the input signal period. 

The basic data collection, conditioning and recording system is also 
shown in Fig. 43. Signals from the single strain gage (per channel) mounted 
on the test piece were conditioned by a strain gage translator. This module 
provided the potentiometric bridge balance resistor, bridge excitation voltage 
and amplification of the normalized strain gage signals to the required sig- 
nal level for recording. By setting the correct K-factor on the input dial 
a normalized signal was obtained equal to 12.5 mv rms per 1000 psi (6.895 MPa) 
peak stress. The factor is given by: 
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0.009065X Ex(RG + RB + RL)^ 
GFx RGxRB 


where RG = Strain gage resistance, ohms 

RB = Balance resistance (=500 ohms) 

RL = Lead wire resistance (ohms) 

E = Youngs Modulus of specimen material (GPa) 

GF = Gage Factor 

The calibrated strain signal level was read directly on a digital voltmeter. 
This was possible since most testing required resonance dwell techniques. From 
the translator, the signal was fed through a Krohnhite variable filter /amplif ier 
and then to a B & K level recorder operating in the peak detection logarithmic 
level mode for permanent chart recording. Alternate recording systems were 
provided by the oscilloscope and the log conver ter/X-Y plotter. Calibration 
of the recorders was performed by recording a signal of known millivolt level, 
as determined by the voltmeter, and using the 20 and 40 dB amplification switches 
on the Krohnhite amplifier to cover the range required. 

Signal conditioning for the response accelerometers was provided either 
by the second channel charge amplifier on the vibration meter or by integral 
electronics and a separate conditioning unit depending on the type of acceler- 
ometer used. Acceleration levels were read off the vibration meter and recorded 
on the level recorder. 


B.3 Data Reduction 

Two types of data plotting were used in the test program, namely, response 
peak logarithmic level decay with time at constant frequency and response 
peak logarithmic level variation with time for time varying input frequency 
(logarithmic sweep about a resonance). From these presentations it was required 
to determine the modal loss factor and frequency data. 

Loss Factor From Decay Curves 


Assuming an equivalent linear viscous damping in a system, the exponential 
envelope of the peaks of a decaying response is given by 


X = XqS 


-7n7t/Tn 


(B.l) 


where x is the amplitude at time t, Xq is the initial or reference amplitude, 
n is the loss factor, and is the period of vibration. The level recorder 
gives the change in amplitude in terms of decibels referenced to Xq. Denoting 
the record level at time t by L^, 
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(B.2) 


Lx = 20 logio(x/xo) 

= “207T (logjo©) T7t/fn = -27.3i7t/Tn 


and 


dLx 

dt 


-27.377/Tn 


(B.3) 


Equation (B.3) can be rewritten as 


dLx _ dLx 
dt " dt 


dt 


= S.P = -27.3T7/Tn 


(B.4) 


where S is the instantaneous slope (dB/mm) of a curve on the chart, £ is the 
length (inm) on the chart and P is the paper speed (mm/sec) of the recorder. 
The loss factor can be calculated by using Eq. (B.4) as 


17 =SPTr^/27 3 


(B.5) 


In practice, the decays were not true straight lines, i.e., q varied with 
strain level and initial transients caused some distortion at the start of the 
decay. An average loss factor was then determined from the slope between 
amplitude levels of 1 and 4 dB down from the initial peak level. Thus, the 
critical region for damping, i.e., the peak to half power level, was spanned. 

Loss Factor From Frequency Response Plots 

The frequency sweep is of the form 


f(t) = fo e“’ 


(B.6) 


where fQ 

f(t) 


is the initial or reference frequency 

is the frequency at time t, and 

is a constant for a given sweep rate. 


89 



Then 


Let At be the time taken for an octave change in frequency, 


f(t-HAt) 

f(t) 


aSt 

e 


From which (B.7) 

a = Cn2/At 

Hence, given the sweep rate (octaves/second), chart paper speed P (mm/sec) and 
a reference frequency f^, the frequency during the sweep can be determined 
from 


f(t) = fo6 


aAC/P 


(B.8) 


where Ai is the measured length (mm) on the chart between fo and f(t). 
The loss factor is given by 


V = 


M. 

tn 


(B.9) 


where Af is the half power bandwidth about the resonant frequency, fn. If 
Af is swept in time At and assuming the response curve is symmetrical about 
the peak, then 


a At/2 , . 

A f = 2(fp 6 

. , a At/2 . 

= 2fn(e -I) 

= fpCiAt for small aAt 


An estimate of loss factor is therefore given by using (B.9) 


7) = a At 


(B.ll) 


Now the measured chart length AJl, corresponding to the change in time At is 
given by 
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Ae = R At 


(B.12) 


where P is the chart paper speed. 

Therefore, the loss factor can be determined using equations (B.ll) and 
(B.12) as 


a AC 

17 =— (B.13) 

In practice, the slowest sweep rate available was chosen, namely 2.2598 x 10”^ 
octaves/sec, giving a value for a of 0.00015664. The loss factors were 
determined by measuring the distance between points 3 dB down from the maxi- 
mum response and taking the resonant frequency as the frequency midway between 
these points. 
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APPENDIX C: DAMPING IN AN ASSEMBLY OF BLADES 


The present program dealt exclusively with single blades and blade-like 
components. An attempt was made to extend the analyses developed, and results 
obtained in the present program to predict damping levels in an assembly. In 
this Chapter, procedures (based on single blade analyses and results) are de- 
veloped for determining the damping levels in an assembly. At the minimum, 
these procedures should yield results which are good approximations to the actual 
levels of damping in an assembly of blades. 

The following assumptions are made. 

1. The assembly is vibrating in a single mode characterized by a nodal 
diameter pattern and a constant interblade phase angle. 

2. The following modal quantities associated with the mode of vibration 

are known: frequency, modal mass, modal stiffness, and modal stress distribution 

of the reference blade (defined as the blade having the maximum amplitude) . 

3. All modal quantities are normalized to unit tip displacement of the 
reference blade. 


7.1 Material Damping 


The procedure for determining the assembly material damping loss factor is 
analogous to the one used in Section 3.1 for a single blade. The loss factor 
can be written as 



N* 

Id 

i=i 


N* 

Z u 


j=i 


(C.l) 


where (*) refers to the assembly, j is the blade number, D the energy dissipated, 
U the strain energy, and N* the number of blades in the assembly. The energy 
dissipated is a function of blade stress and damping constants. The strain energy 
is a function of blade stress and modulus of elasticity. The stress at a point 
on the j can be written as 

O'. = o-p COS j /G*=27rp/N* , (C.2) 
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where o is the stress at the corresponding point on the reference blade, p is 
the number of nodal diameters, and 3* is the interblade phase angle. Equations 
(C.2) and (3.4) lead to 


Dj = Dr(C 0 S j 0* , U| = Ur (cos 

where n is the damping constant. 


(C.3) 


Substitution of Eq. (C.3) into (C.l) and introduction of 


Z cos^ 

j=l 

Z cos^ ]j3 


Cn ~ 


l_ Pr 

'^R " 27T Ur * 


(C.4) 


leads to 




(C.5) 


where q is the material damping loss factor for the reference blade, computed by 
using tSe blade modal stress distribution for the assembly mode of vibration and 
the analysis of Section 3.1. Clearly, for n = 2, C = 1 so that the loss factors 
for a reference blade and an assembly are identical. 


7 . 2 Friction Damping 

The friction damping effects considered in this section are those due to 
slipping at blade roots and rubbing at shroud interfaces. An assembly mode of 
vibration which is characterized by modal quantities can be thought of as a single 
degree of freedom, acted on by friction forces. Thus, the possibility of using 
the analyses developed for single blades, with some modifications, to predict the 
effects of friction damping in an assembly of blades is explored here. In part- 
icular, the idea is to use the computer programs ROTDMP and SHDDMP developed for 
single blades, with suitable modifications to the input information to reflect 
the dynamics of the vibrating assembly. Brief descriptions of the necessary 
modifications to the input information are as follows. 

(1) Modal Information: 

The modal parameters; frequency f’’^, modal stiffness k*, and modal mass m* 
are those that characterize the assembly mode of vibration. 

(2) Baseline Condition: 

The baseline condition is assumed to be condition of no friction damping. 
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The baseline stress a is the blade modal stress at some point on the reference 

blade. The baseline damping loss factor p * is assumed to be known, 

o 

(3) Base Excitation Level in g-units , n 

§ 

Typically, this information is unavailable for an assembly of blades. 
However, Ug can be estimated provided some additional information pertaining to 
assembly is given as follows: 

(a) Given , the excitation force amplitude, 



(2TTf*f 

g 


(C.6) 


where g is the acceleration due to gravity 

(b) Given x *, baseline tip deflection of 
o 


ng 


[Z-rrff 
Cm Q 


(9.81 m/sec^) 
reference blade. 


(C.7) 


where G is a constant depending on the mode of vibration. For example, 

for a uniform cantilever beam C = 1.566 for the first mode (Ref. 39). 

However, in the case of an assembly C may have to be determined from 

m 

tests . 

(4) Interface normal load, coefficient of friction: 

In single blade analysis the energy dissipated due to friction at an interface 
is computed as. 


D = 4 /J. N Xg » 


(C.8) 


where N is the normal load, y the coefficient of friction, and Xg the slip 
amplitude. In an assembly, all the blades are not equally active due to the 
existence of a nodal diameter pattern. The total energy dissipated in the assem- 
bly can be estimated as 

* N* 

= 4yU N Xg Z 

j=l 


cos J p 


(C.9) 


Thus, a convenient way of accounting for the dynamics of the assembly is to define 
an equivalent coefficient of friction. 
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(C.IO) 



= Z cos j /3* . 
i=i 


(5) Mass Ratio, 3 (Root Damping) 

For a single blade the mass ratio is defined as 


/S = (Iq/L^ )/m| 


(C.ll) 


where I is the blade mass moment of inertia about the root, L the blade length, 
o 

and m^ the modal mass. For an assembly, the total mass moment of inertia about 
the disk circumference can be written as 



so that, an equivalent mass ratio for the assembly is, 

/ 8 *= . 


(C.12) 


(C.13) 


(6) Stiffness ratio 3, (Shroud Damping) 

For a single blade the stiffness is defined as 

P = R/k , 


(C.14) 


A 

where k is the static stiffness of the blade below the shrouds. In an assembly, 
this static stiffness is not the same as that for a single blade due to the in- 
fluence of neighboring blades. However, the appropriate value of the static 
stiffness, can be obtained from static tests on the assembly, so that an equivalent 
stiffness ratio can be defined as 


/8 


* 


A* * 

k /k 


(C.15) 
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TABLE I 


TITANIM CANTILEVER BEAM TEST RESULTS 


Input Level^ 

Mode g pk 

1 0.010 

0.012 

0.014 

0.015 

0.0165 

0.018 

0.020 

0.022 

0.023 

0.025 

2 0.05 

0.08 

0.10 

0.15 

0.16 

0.20 

0.20 

0.22 

0.24 

0.25 

0.26 

0.28 

0.30 

0.32 

0.34 

0.35 

0.36 

0.40 

3 0.10 

0.20 
0.20 
0.30 
0.40 
0.50 
0.55 
0.60 
0.70 
0.80 
1.00 
1 20 

NOTES: 


Period 

Response Level 

sec 

microstrain 

0.0146220 

511 

0.0146216 

662 

0.0146158 

702 

0.0146160 

822 

0.0146200 

933 

0.0146196 

1076 

0.0146190 

1178 

0.0146187 

1276 

0.0146206 

1293 

0.0146186 

1320 

0.0023390 

138 

0.0023392 

222 

0.0023390 

262 

0.0023391 

409 

0.0023395 

427 

0.0023392 

498 

0.0023397 

520 

0.0023399 

591 

0.0023398 

609 

0.0023395 

609 

0.0023401 

680 

0.0023401 

716 

0.0023391 

733 

0.0023401 

778 

0.0023402 

827 

0.0023402 

911 

0.0023401 

867 

0.0023405 

1093 

0.00083618 

71 

0.00083598 

147 

0.00083622 

147 

0.00083618 

222 

0.00083610 

293 

0.00083636 

369 

0.00083620 

400 

0.00083630 

440 

0.00083651 

507 

0.00083649 

564 

0.00083675 

680 

0.00083709 

791 


ISinusoidal vibration applied at root of cantilever 
^Axial strain measured at 8% span 
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TABLE II 


SUMMARY OF TWISTED PLATE TEST RESULTS 

Vibratory 

Nominal Centr. Stress^ to Static 

Rotational Force Mode Freq., Range -Stress Loss Factor 

Speed, rpm N No. _ Hz MPa (pk) Ratio Range Average 

Centrifugal Test 



(1 

100.7 

16.2-27.6 


- 

.0014 _ 

.0017 

.0016 

0 

0 {2 

549 

27.0-39.4 


- 

.00040 - 

.00071 

.00059 

(in air) 

(3 

677 

35.3 


- 

.00030 - 

.00031 

.00031 


(1 

101 

27.6 


— 

.00049 - 

.00051 

.00050 

0 

0 ^2 

550 

27.6 


- 

.00040 - 

.00042 

.00040 

(in vac) 

(3 

677 

39.0 


- 

.00015 - 

.00017 

.00016 


(1 

103 

13.8-27.6 

8.3 

16.6 

.00050 _ 

.00071 

.00060 

500 

211 2 

553 

27.6-33.0 

16.6 

19.9 

.00055 - 

.00065 

.00060 


(3 

678 

30.0-39.6 

18.1 

23.9 

.00020 - 

.00022 

.00020 


(1 

109 

13.8-27.6 

2.1 

4.2 

.00055 - 

.00102 

.00080 

1000 

845 <2 

559 

11.0-46.0 

1.7 

6.9 

.00060 - 

.00079 

.00070 


(3 

677 

13.0-24.8 

2.0 

3.7 

.00019 - 

.00022 

.00021 


(1 

113 

27.6-30.9 

2.9 

3.2 

.00093 - 

.00104 

.00097 


1200 1220 ^2 

(3 



(1 

117 

9.9 

0.76 

.00080 _ 

.00104 

.00093 

1400 

1650 {2 

566 

9.9-19.9 

0.76 

1.5 .00067 _ 

.00085 

.00074 


(3 

678 

7.2-13.2 

0.55 

1.0 .00019 - 

.00025 

.00022 


(1 

119 

8.8-11.6 

0.59 

0.78.00085 - 

.00111 

.00099 

1500 

1900 <2 

568 

9.9-12.7 

0.66 

0.85.00040 - 

.00060 

.00051 


(3 

678 

9.9-22.2 

0.66 

1.5 .00018 - 

.00023 

.00020 


2 

Bench Test 



(1 

101 

9.3-86 

- 

.00056 _ 

.00088 

.00070 

(in air) 

0 2 

548 

25-97 

- 

.00076 - 

.00107 

.00097 


(3 

679 

5.8-30 

- 

.00013 ^ 

.00014 

.00014 


(1 

101 

66-104 

— 

.00049 _ 

.00054 

.00052 

(in vac) 

0 < 2 

549 

21-95 

- 

.00027 

.00036 

.00032 


(3 

680 

7.9-41 

- 

.000092 - 

.0000104 .000097 

NOTES ; 1 

Stress value 

assuming E - 

128 GPa 






2 Steel shim removed from tip 
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TABLE III 


Pressure 

Ratio^ 

P/Ps 


0.02 


0.26 


1.0 


NOTES : 


SUMMARY OF DATA FOR FAN BLADE MATERIAL DAMPING TESTS 


Input 

Response 

Root 

ASMT 


Acceleration 

Frequency 

Stress 

Stress 

Loss 

. R 

Mode Hz 

MPa 

MPa 

Factor 



1(1F) 

4 




ND 

77.88 

56 

28 

.00030' 

ND 


77.89 

110 

52 

.00085 

0.19 

2(2F) 

250.38 

47 

32 

.00024 

0.56 


250.39 

145 

98 

.00027 

0.83 


250.37 

207 

141 

.00027 

0.6 

3(1T) 

483.09 

36 

47 

.00011 

1.0 


482.91 

58 

76 

.00012 

2.8 


481.39 

126 

165 

.00014 

6.0 


480.35 

159 

207 

.00015 

ND 

1(1F) 

77.80 

56 

ND 

.0013 

ND 


77.77 

110 

ND 

.0023 

0.38 

2(2F) 

250.29 

68 

ND 

.00032 

0.85 


250.29 

128 

ND 

.00049 

1.7 


250.29 

200 

ND 

.00062 

2.5 

3(1T) 

482.46 

65 

ND 

.00029 

7.0 


481.19 

123 

ND 

.00048 

ND 

1(1F) 

77.57 

56 

ND 

.0040 

ND 


77.44 

110 

ND 

.0062 

0.75 

2(2F) 

250.09 

62 

ND 

.00067 

2.0 


249.96 

133 

ND 

.00109 

4.0 


249.96 

196 

ND 

.00130 

5.7 

3(1T) 

481.72 

53 

70 

.00087 

7.0 


481.74 

63 

82 

.00095 

10.0 


481.58 

72 

95 

.00102 


1. Pg = 760 Torr 

2. Loss factors for mode 1 derived from frequency sweep data 

3. Loss factors for modes 2 and 3 derived from decay data 

4. Frequencies estimated from Figure 20 
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TABLE IV 


SUIIMARY OF ROOT DAMPING TEST DATA 


Total 

Excitation 

Root 

Response 

Bolt Load 

Level 

Stress 

Frequency 

N 

g(pk) 

MPa (pk) 

Hz 

576 

0.2 

5.4 

68.311 


0.2 

5.5 

68.353 


0.5 

9.9 

66.445 


0.5 

9.9 

66.454 


1.0 

11.0 

64.135 


1.0 

11.0 

64.144 

1152 

0.2 

6.1 

72.396 


0.2 

6.1 

72.369 


0.2 

6.1 

72.369 


0.5 

12.1 

69.989 


0.5 

12.4 

70.013 


1.0 

17.4 

68.871 


1.0 

17.4 

68.890 


1.0 

17.4 

68.899 

1727 

0.2 

11.5 

73.910 


0.2 

11.2 

73.915 


0.2 

11.2 

73.893 


0.5 

14.8 

72.516 


0.5 

15.2 

72.490 


0.5 

15.4 

72.543 


1.0 

21.2 

71.367 


1.0 

21.2 

71.408 

2303 

0.2 

18.0 

74.460 


0.2 

17.8 

74.555 


0.2 

18.1 

74.560 


0.5 

23.9 

73.719 


0.5 

24.4 

73.735 


0.5 

24.4 

73.719 


1.0 

28.8 

72.690 


1.0 

28.8 

72.690 


1.0 

28.8 

72.690 

2879 

0.2 

21.9 

74.794 


0.2 

21.9 

74.800 


0.5 

31.0 

74.173 


0.5 

30.9 

74.184 


1.0 

38.4 

73.341 


1.0 

38.9 

73.341 


1.0 

38.8 

73.346 


Loss 

Factor 

n 

.0374 

.0374 

.0501 

.0501 

.0846 

.0815 

.0259 

.0234 

.0259 

.0436 

.0436 

.0515 

.0550 

.0469 

.0129 

.0132 

.0135 

.0289 

.0285 

.0273 

.0446 

.0433 

.0089 

.0088 

.0088 

.0149 

.0135 

.0147 

.0263 

.0265 

.0284 

.0073 

.0073 

.0113 

.0113 

.0207 

.0208 

.0198 
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TABLE V 


SUMMARY OF SHROUDED BLADE RESPONSE TEST DATA 


Shroud 

Input 

Excit . 

ASMT 

Loss 

Shroud 0.0 

.P.2 

Shroud 

I.P.^ 

Normal 

Accel . 

Period 

Stress 

Factor 

Acceleration g 

Acceleration, g 

Load N 

g(pk) 

ms 

MPa(pk) 

ri 

Unfilt^ 

Filt. 

Unfilt 

Filt, 




First 

Flap Mode 

Results 




150 

0.5 

3.566 

41 

.0034 

2.5 

2.5 

1.8 

1.2 


1.0 

3.574 

54 

.0038 

2.9 

2.9 

2.7 

1.8 


2.0 

3.589 

59 

.0056 

3.6 

3.6 




3.0 

3.589 

60 

.0054 






4.0 

3.579 

63 

.0053 





300 

0.5 

3.563 

61 

.0024 

3.8 

3.8 

1.7 

1.6 


1.0 

3.566 

103 

.0023 

7.2 

7.0 

2.1 

1.7 


2.0 

3.574 

124 

.0036 

8.8 

8.8 

4.0 

3.2 


3.0 

3.575 

128 

.0044 

9.6 

9.1 

5.2 

2.6 


4.0 

3.580 

134 

.0036 

10.0 

7.6 

9.0 

3.6 

500 

0.25 

3.558 

39 

.0031 






0.5 

3.551 

56 

.0026 

4.0 

4.0 




0.5 

3.558 

60 

.0024 



1.5 

1.4 


1.0 

3.556 

102 

.0027 

6.8 

6.8 




1.0 

3.559 

124 

.0022 






2.0 

3.560 

154 

.0034 

12 .0 

11 .0 




2.0 

3.564 

163 

.0031 



4.7 

3.1 


3.0 

3.562 

178 

.0031 

13.0 

11.0 




3.1 

3.566 

182 

.0034 



8.2 



4.0 

3.565 

186 

.0031 

14.0 

14 .0 




4.0 

3.565 

195 

.0032 








First Torsion Mode 

Results 




150 

0.5 

0.8932 

5.5 

.00052 

2.6 

2.1 

2.0 

1.8 


1.0 

0.8940 

9.2 


3.8 

3.6 

3.0 

2.7 


2.0 

0.8944 

12 .0 

.00080 

5.6 

5.0 

1.3 

0.9 

300 

0.5 

0.8932 

6.1 

.00058 

3.4 

3.4 

0.9 

0.8 


1.0 

0.8935 

11 .0 

.00055 

5.8 

5.2 

2.4 

1.2 


2.0 

0.8933 

15 .0 

.00075 

4.0 

3.6 

3.0 

1.3 


3.0 

0.8940 

20 .0 

.00082 

14 .0 

9.9 

9.9 

3.5 

500 

0.5 

0.8910 

7.7 

.00040 

4.4 

3.8 

1.5 

1.4 


1.0 

0.8912 

13 .0 

.00045 

6.2 

6.2 

1.9 

1.7 


2..0 

0.8918 

17 .0 

.00059 

8.4 

8.4 

3.8 

3.0 


3.0 

0.8924 

23 .0 

.00074 


11 .0 

3.2 

3.1 


4.0 

0.8933 

28 .0 

.00084 

14 .0 

13 .0 

3.2 

1.8 

Notes: 

1. Slop 

e of decay 

' taken 

between the 

1 and 4dB 

points 




2. 0.0. 

P = Out of 

■ plane 

(normal to 

shroud plane) 




3. I.P 

= In plane 

i (normal to shroud plane) 





4. Low pass filter 
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TABLE VI 


SUMMARY OF FRICTION TEST RESULTS 


Input 

Input 

Relative 

Maximum 

Normal 

Work Done 

Equivalent 

Accel 

Freq. 

Disp. 

Velocity 

Load 

Per Cycle 

Friction 

g 

Hz 

mm(pk-pk) 

mm/ sec(pk) 

N 

mJ 

Coefficient 

1.7 

80 

0.132 

33.2 

271 

14.6 

0.407 

1.7 

80 

0.132 

33.2 

494 

22.3 

0.341 

1.7 

80 

0.132 

33.2 

672 

28.1 

0.317 

3.0 

140 

0.076 

33.5 

271 

10.7 

0.517 

3.0 

140 

0.076 

33.5 

494 

15.7 

0.417 

3.0 

140 

0.076 

33.5 

672 

18.9 

0.370 

5.0 

140 

0.127 

55.9 

271 

17.5 

0.509 

5.0 

140 

0.127 

55.9 

494 

27.2 

0.434 

5.0 

140 

0.127 

55.9 

672 

33.7 

0.395 

10.0 

280 

0.064 

55.9 

271 

7.7 

0.445 

10.0 

280 

0.064 

55.9 

494 

11.5 

0.378 

10.0 

280 

0.064 

55.9 

672 

14.0 

0.329 

20.0 

280 

0.127 

112 

271 

16.6 

0.482 

20.0 

280 

0.127 

112 

494 

24.4 

0.390 

20.0 

280 

0.127 

112 

672 

30.4 

0.357 

20.0 

420 

0.056 

73.7 

271 

7.4 

0.487 

20.0 

420 

0.056 

73.7 

494 

10.6 

0.382 

20.0 

420 

0.056 

73.7 

672 

13.1 

0.345 
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TABLE VII 


SUMMARY 

OF TURBINE BLADE 

RESPONSE TEST 

RESULTS 

Input 

Input 

Shank 

Root 

Level 

Frequency 

Stress 

Stress 

g(pk) 

Hz 

MPa(pk) 

MPa (pk) 

0.06 

515.498 

34 

25 

0.10 

515.328 

53 

42 

0.20 

515.291 

97 

79 

0.30 

515.286 

133 

108 

0.40 

515.129 

139 

no 

0.50 

515.132 

181 

146 

0.60 

515.137 

197 

158 

0.70 

515.137 

209 

169 

0.79 

515.145 

220 

179 

0.9 

515.002 

236 

191 
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TABLE VIII 


SUMMARY OF TURBINE BLADE PLATFORM DAMPING TEST RESULTS ("SMALL DAMPER") 


Normal 

Input 

Input 

Shank 

Root 

Normal 

Input 

Input 

Shank 

Root 

Load 

Level 

Frequency 

Stress 

Stress 

Load 

Level 

Frequency 

Stress 

Stress 

N 

g(pk) 

Hz 

MPa(pk) 

MPa(pk) 

N 

g(pT^) 

Hz 

MPa(pk) 

MPa(pk) 

1.65 

0.29 

515.892 

41 

30 

26.4 

3.9 

516.025 

44 

37 


0.40 

515.573 

109 

88 


4.1 

515.868 

259 

204 


0.50 

515.591 

136 

111 







0.60 

515.467 

149 

120 

33 

3.0 

515.198 

12 

- 


0.70 

515.331 

164 

135 


4.0 

516.377 

26 

22 


0.80 

515.374 

178 

138 


5.0 

516.649 

88 

67 


1.0 

515.294 

193 

157 


6.0 

516.646 

>300 

>300 

3.3 

0.55 

515.969 

99 

79 

88.9 

1 

562.212 

18 

19 


0.70 

515.498 

153 

123 


2 

561.325 

29 

30 


0.85 

515.464 

168 

137 


3 

558.438 

39 

41 


1.0 

515.315 

194 

157 


4 

560.601 

38 

40 







5 

561.558 

38 

40 

4.9 

0.70 

515.921 

72 

- 


6 

563.304 

38 

40 


0.77 

516.396 

66 

47 


7 

566.640 

35 

36 


0.85 

515.871 

135 

112 


8 

568.971 

33 

34 


1.0 

515.565 

175 

152 


8 

515.039 

47 

40 







9 

515.464 

55 

47 

6 . 6 

1.1 

515.783 

129 

101 


9 

577.247 

26 

31 


1.5 

515.778 

200 

160 


10 

575.530 

29 

33 







10 

515.198 

56 

50 

13.2 

1.8 

515.823 

7 

- 







2.0 

515.879 

215 

182 






19.8 

2.6 

515.464 

7 

7 

178 

2 

565.576 

20 

25 


2.8 

515.464 

12 

10 


4 

559.059 

45 

52 


3.0 

515.703 

210 

182 


4 

565.323 

65 

73 







6 

558.090 

77 

85 







8 

557.302 

85 

93 







10 

554.416 

89 

99 



104 


TABLE IX 


SUMMARY OF TURBINE BLADE PLATFORM DAMPING TEST RESULTS ("LARGE” DAMPER) 


Normal 

Load 

N 

Input 

Level 

g(pk) 

Input 

Frequency 

Hz 

Shank 

Stress 

MPa(pk) 

Root 

Stress 

MPa(pk) 

Normal 

Load 

N 

Input 

Level 

g(pk) 

Input 

Frequency 

Hz 

Shank 

Stress 

MPa(pk) 

Root 

Stress 

MPa(pk) 

1.65 

-0.255 

517.786 

39 

28 

19.8 

2.1 

516.486 

20 

ND 


0.28 

516.028 

41 

30 


2.2 

516.134 

160 




-0.30 

515.903 

73 

60 


2.3 

516.129 

215 




0.315 

515.730 

79 

63 








-0.34 

515.735 

83 

69 

26.4 

2.8 

ND 

11 




0.36 

515.730 

96 

77 


3.0 

516.145 

18 




-0.40 

515.663 

108 

88 


3.2 

515.951 

243 




-0.45 

515.695 

119 

94 








-0.51 

515.453 

140 

114 








0.55 

515.419 

155 

125 








0.60 

515.382 

144 

118 








0.65 

515.347 

154 

118 

88.9 

2 

557.072 

17 

ND 


0.70 

515.305 

154 

123 


3 

540.003 

22 

ND 


0.80 

515.315 

173 

137 


4 

524.084 

24 

ND 








5 

510.618 

28 

21 

3.3 

0.45 

516.630 

42 

ND 


6 

511.546 

35 

ND 


0.55 

516.332 

73 

ND 


7 

511.339 

36 

ND 


0.70 

515.573 

143 

116 


8 

507.261 

47 

38 








9 

515.265 

>300 

ND 

6.6 

0.85 

515.959 

113 

ND 


10 

516.156 

>300 

ND 


0.90 

515.919 

126 










1.0 

515.969 

147 



178 

2 

584.669 

23 

28 








3 

581.727 

40 

ND 

9.9 

0.9 

516.140 

79 




4 

580.649 

47 

ND 


1.1 

516.165 

140 




5 

574.511 

44 

ND 


1.3 

516.134 

164 




6 

574.511 

52 

61 








7 

574.260 

52 

59 

13.2 

1.5 

516.097 

93-124 




8 

ND 

50 

57 


1.9 

516.129 

195 




10 

574.548 

56 

61 



TABLE B-1 


TEST COMPONENTS AND OTHER HARDWARE DESIGNED 
AND FABRICATED FOR THE TEST PROGRAM 


Item 

No. 


Item Description 

Material Constants Determination 

1. Test Piece (Uniform Cantilever Beam) 

2. Adapter Clamp 

3. Steel Root Block 

4. Vacuum Chamber 


Design and Fabrication Refs. Section 
Drawing Nos/Engineering Order Nos . Ref . 


EOI513005 
EOI513003&4 
As For Item V-2 
As For Item III-2 


3.3.1 


II Twisted Plate-Centrifuge & Bench Tests 

1. Test Piece (Twisted Plate) 

2. Support Assembly 

3. Driver Magnet Assembly 

4. Adapter Clamp 

5. Steel Root Block 

6. Vacuum Chamber 


3.3.2 


DRG 1794-1 Parts List 48300 
+ EOI513000 

Unlisted DRG DTD 7-31-79 
EOI513007 
As For Item V-2 
As For Item III-2 


III Fan Blade-Evacuated & Thermal Tests 

1. Test Piece (Blade Welded in Titanium 

- Root Block) 

2. Vacuum Chamber 

3. Thermal Lid 


(Blade Furnished by P&WA) 
DRG1794-2 Parts List 48466 


3.3.3 


IV Shrouded Fan Blade Tests 

1. Test Piece (Fan Blade Welded at Root) 

2. Load Platens 

3. Vacuum Chamber 

4. Load Mechanism Assembly 

V Root Damping Tests 

1. Test Piece (Fan Blade) 

2. Steel Root Block 

3. Root Attachment Assembly 

4. Load Measurement Bolts 


5.2.1 


(Blade Furnished by P&WA) 
DRG1794-2 Parts List 48466 


(Blade Furnished by P&WA) 
DRG1794-2 Parts List 48466 
+ EOI513009 


4.2.1 


VI Platform Damping Tests 

1. Test Piece (Turbine Blade) 

2. Fixture-Root Clamp and Base Plate 

3. Loading Arm Includes Damper (2 Made) 

VII Composite Beam Damping 

1. Test Piece (Flat Strlp-4 Specimens) 

2. Adapter Clamp 

3. Steel Root Block 

4. Vacuum Chamber 


6 . 2.1 


(Blade Furnished by P&WA) 
DRG1794-3 Parts List 48595 


3.3.4 

Furnished by UTRC S/N 3491A/B&3492A/B 
Unlisted DRG 
As For Item V-2 
As For Item III-2 


VIII Blade Coating Investigation 

1. Test Piece (Turbine Blade) 

2. Fixture 


3.3.5 


As For Item VI-1 Stripped & Coated 
As For Item VI-2 


IX Friction Evaluation Tests 

1. Test Piece (1 Pair Load Platens) 

2. Fixture Assembly 

3. Load Measurement Bolts 


5.3.1 


As For Item IV-2 
E0I513017 thru 021 
As For Item V-4 


105 



TABLE B-2 


EQUIPMENT USED IN TEST PROGRAM 
Model or 


Item No. 

Equipment 

Type No. 

Serial No. 

Manufacturer 


Environmental, Load 

and Input Excitation 

Equipment 


1 

Electro-Dynamic Vibrator (5300 N) 

724 

128 

L.D.S 

2 

Power Amplifier 2 kW (used with Item 1) 

MPA 2 

139 

L.D.S 

3 

Power Amplifier 400 W 

2250 MB 


MB Electronics 

4 

Magnetic Driver 75 W 

- 

- 

Built In-house 

5 

Centrifuge (1.73 m dia, 3000 rpm) 

- 

_ 

UTRC Rig 

6 

Vacuum Pump 

3164 A 

- 

Kinney 

7 

Load Bolts (1500 N capacity) 

- 

- 

Built In-house 

8 

Automatic Vibration Exciter Control 

1025 

199899 

B & K 

9 

Vibration Meter 

DVA 

270 

L.D.S 

10 

Dynamic Analyser 50 Hz BW Tracking Filter 

SDIOIB 

84 

Spectral Dynamics 

11 

Variable Filter/Amplifier (40 dB) 

3322 

139 

Krohnhite 

12 

Variac (2) 

W5 


G. R. 

13 

TC Readout Meter 

2165 A 

0960046 

Fluke 

14 

Vacuum Gage 

FA160 

YY10223 

Wallace & Tiernan 

15 

Portable Strain Indicator 

P350 

2486 

Budd 

16 

Strain Gage Bridge Amplifier 

BAM IB 

14194 

Vishay Instr. 

17 

Counter/Timer 

1953 A 

2195003 

Fluke 

18 

Counter/Timer 

5532A 

209-00170 

H.P. 

19 

Accelerometers 

101 

105/106 

V.S.S. 

20 

Strain Gages (on load bolts - Item 7) 

EA-06-050-AH-120 


MM 

21 

Thermocouples Chromel/Alumel 

- 

- 

Built In-house 


Response Data Collection, 

Conditioning and Recording Equipment 


22 

Accelerometers; Sensitivity 41.8/46.5 pcmb/g 

101 

105/106 

V.S.S. 

23 

Sensitivity 10 mv/g 

MIOOO 

6008 

Vibrametrics , Inc. 

24 

Sensitivity 5 mv/g 

2250 

N195 

Endevco 

25 

Strain Gages: Miniature Foil Type 

EA-05-062-AK-120 


MM 

26 

Miniature Foil Type 

EA-06-062-AK-120 


MM 

27 

Mid Temperature Type 

FSM-12-35-S6 


BLH 

28 

Vibration Meter 

As for Item 9 


L.D.S. 

29 

Signal Conditioner (used with Item 23) 

LP15-3 

393 

Vibrametrics, Inc. 

30 

Signal Conditioner (used with Item 24) 

4416 

AB27 

Endevco 

31 

Strain Gage Translators 

ED1127-8 

128/133 

Built In-house 

32 

Strain Gage Bridge Amplifier 

As for item 16 



33 

Variable Amplifier Filter 

As for Item 11 



34 

Log Converter 

7562A 

1211A01693 

H.P. 

35 

Level Recorder 

2305 

240896 

B & K 

36 

Digital Voltmeters 

8000A 

0645437 

Fluke 

37 

Counter /timer 

As for items 17 

& 18 


38 

X-Y Recorder 

7035B 

1722A13849 

H.P. 

39 

Oscilloscope with camera 

533 

001791 

Tektronix 

40 

Time Base Unit (used with Item 38) 

17108A 

938A01100 

H. P. 

41 

Digital Multimeter 

7004A 

2878 

Systron Donner 
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ORIGINAL PAGE IS 79-09 
OF POOR QUALITY 






ROTATION 

DIRECTION 


Figure 3. Twisted Plate Installation for Centrifugal Test 
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Figure 5. Hysteresis Loop for Material Damping 




CASE 1 




I / 

z=o — — 


CASE 2 

Figure 6b. Stress Distribution Across Element Thickness 
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EXCITATION 


STRAIN GAGE 



7a. GENERAL ARRANGEMENT OF TEST SET-UP 


VACUUM VIBRATION FIXTURE 



TOP VIEW END VIEW 


7b. TITANIUM (8-1-1 ) SINGLE CANTILEVER BEAM TEST SPECIMEN 

Figure 7. Material Damping Constant Measurement Test Hardware 
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Figure 9. Twisted Plate Mounted in Bench Test Fixture 
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Figure 10. NASTRAN Finite Element Model for the Twisted Plate 
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Figure 11. Biock Diagram of the Driving and Recording Systems used to Obtain Vibration 
Decay Data for the Twisted Plate on the Centrifuge 


79 - 09 - 64-1 


116 





LOSS FACTOR r] 




ROTOR SPEED, rpm 


Figure 12. Variation of Modai Damping (Loss Factor) with Rotationai Speed for 
Titanium Twisted Plate in Partial Vacuum (5 torr) 
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RUN 001 14V 


MODE 1 (IF) 

FREQUENCY ==101,1 Hz 
MAX STRESS = 27,6 MPa 
LOSS FACTOR = 0,00049 


WRITING SPEED 50 mm/sec 
PAPER SPEED 3 mm/sec 



RUN 002 17V 


MODE 2 (2F) 
FREQUENCY = 550,3 Hz 
MAX STRESS = 27,6 MPa 
LOSS FACTOR = 0.00042 


WRITING SPEED 200 mm/sec 
PAPER SPEED 10 mm/sec 


MODE 3 (IT) 
FREQUENCY = 677.2 Hz 
MAX STRESS = 27.6 MPa 
LOSS FACTOR = 0.0001 6 


WRITING SPEED 200 mm/sec 



PAPER SPEED 10 mm/sec 


Figure 13. Decay Characteristics of Twisted Piate Vibrating in its First 
Three Modes at Zero Speed in Vacuum 
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VIBRATION PEAK LEVEL 



MODE 1 (IF) 
FREQUENCY = 109.2 Hz 
MAX STRESS = 27,6 MPa 
LOSS FACTOR = 0.0007 
SPEED = 1000 rpm 


WRITING SPEED 31.5 mm/sec 
PAPER SPEED 3 mm/sec 



MODE 1 (1F) 
FREQUENCY = 109.4 Hz 
MAX STRESS = 13.8 MPa 
LOSS FACTOR = 0.001 
SPEED = 1011 rpm 


r 


cm 

T 



TIME - 


WRITING SPEED 50 mm/sec 
PAPER SPEED 3 mm/sec 


Figure 14. Decay Characteristics of Twisted Piate Vibrating in its First Mode 
at Speed in Vacuum. (Two Samples) 
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VIBRATION PEAK LEVEL (dB) 


MODE 2 (2F) 


: RUN 102 61V 



FREQUENCY = 557.5 Hz 
MAX STRESS = 25.4 MPa 
LOSS FACTOR = 0.00066 
SPEED = 978 rpm 




cm 


~r - - p- -T 1 - - - 


WRITING SPEED 31 .5 mm/sec 
PAPER SPEED 10 mm/sec 


MODE 3 (IT) 


5 mvf RUN 103 82V L 



TIME -* 


Figure 15. Decay Characteristics of Twisted Plate Vibrating in Its Second and 
Third Modes at Speed in Vacuum. (5 torr) 
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Figure 16. Instrumented Fan Blade Showing “Root” Strain Gage 
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Figure 17. Instrumented Fan Blade Showing A.S.M.T. Strain Gage 
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TITANIUM FAN BLADE WELDED AT ROOT- 


HEATING ELEMENTS 


EXCITATION I 
DIRECTION 


ALUMINUM FIXTURE- 


■ TITANIUM BLOCK 


•INSTRUMENTATION PORT 


FIBERGLASS COVER-^ 


INPUT CONTROL 
ACCELEROMETER 


■SHAKER SLIP PLATE 


DOWEL PIN 


II 

)( 

rr^- 


uA- 


nu' 




ROOT THERMOCOUPLE- 


MIDSPAN THERMOCOUPLE- 


TIE DOWN' 
BOLT 


TIP THERMOCOUPLE 


Figure 18. Test Set-Up for Fan Blade Material Damping Investigation at Elevated Temperatures 




NOTES; 

Ps = 760 TORR 

MAXIMUM OF STRESSES MEASURED AT BLADE ROOT (E = 124 GPa) 

FREQUENCY SWEEP RATE = MINIMUM AVAILABLE (0.01 36 OCTAVE PER MINUTE) ALLOWS <5% 

DEVIATION FROM STEADY STATE FOR T? >0.0022 


Figure 19. Variation of Measured Modal Damping with Air 
Pressure for Fan Blade 
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NOTES: 

Ps = 760 TORR 

fs = 77.5795 Hz (MEASURED FREQUENCY AT ATMOSPHERIC PRESSURE, 755 TORR) 

O = MAXIMUM STRESS MEASURED AT BLADE ROOT (E = 124 GPa) 

Figure 20. Variation of Measured First Mode Frequency Ratio with Air Pressure 

for Fan Blade 
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LOSS FACTOR 


0.005 


0.004 


0.003 


0.002 


0.001 


0 



0 10 20 30 40 50 60 


PEAK ROOT STRESS, MPa 


Figure 22. Variation of Loss Factor with Root Stress and Soak Temperature 
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80-220 



Figure 26. Composite Material Cantilevered Beam Vibration Test Set Up 



MEASURED LOSS FACTOR 


KEY; 


KEY 


OOMODE 1, FIRST SET i 
• •MODE 1, SECOND SET f 
+ +MODE 2 


f = 109.2 — 109.8 Hz 
f = 681.2 — 686.3 Hz 


OO MODE 1, FIRST SET > 

• • MODE 1, SECOND SET f ^-99-3 — 99.5 Hz 
+ + MODE 2 ■ f = 620.7 Hz 




(a) B/AI (b) B/AI/Ti 

Figure 27. Composite Material Damping 
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LEGEND: 


(PRINCIPAL STRESS/MAX STRESS) 
NEAR MID-CHORD 
ELEMENT #2 


UPPER SURFACE 

LOWER SURFACE 





(ROOT) SPAN WISE STATION (TIP) 

(CENTER OF ELEMENT) 


Figure 28. Spanwise Distributions of Outer Fiber 

Major Principal Stresses for Twisted Plate 
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ROOT PEAK STRESS 



0 0.05 0.1 0.1 0.2 MODE 2 

INPUT ACCELERATION LEVEL — g 

Figure 29. Twisted Plate Forced Response in Partial Vacuum (10 Torr) 
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Figure 30a. Motions of Blade with Slip at Root 



rrmr mm 


Figure 30b. Analytical Model of Blade with Slip at Root 
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79 - 09-64 


LO 

a^ 


LOADING 

SCREW 


STEEL 

FIXTURE 

BLOCK 



NOTE; FOR TESTING THIS ASSEMBLY IS BOLTED 

INTO THE VACUUM CHAMBER/TEST FIXTURE 
USED FOR FAN BLADE MATERIAL AND SHROUD 
DAMPING TESTS 


DIRECTION 

OF 

EXCITATION 



Figure 31. Test Fixture and Root Loading Device for Fan Blade Root Damping Tests 

I 





Figure 32. Fan Blade/Fixture Test Assembly for Root Damping investigation 

80-69- A 81-3-7-22 
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PORTABLE BRIDGE 
BALANCING AND 
STRAIN INDICATOR UNIT 
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ROOT STRESS (PEAK) — MPa 




EXCITATION LEVEL :=0.2g 
FREQUENCY = 68.31 Hz 
LOSS FACTOR = 0,037 


^-11 MPa 


EXCITATION LEVEL = 1 .Og 

FREQUENCY = 64.1 4 Hz 

.v.: loss factor = o.oas 


■PAPER SPEED lOOmm/sec 


Figure 34. Decays of Blade Response with an Applied Root Loading of 576N. 
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ROOT STRESS (PEAK) — MPa 



EXCITATION LEVEL = 0.2g 
FREQUENCY = 74.79 Hz 
LOSS FACTOR = 0.0072 



TIME 


EXCITATION LEVEL = 1.0g 
FREQUENCY = 73.34 Hz 
LOSS FACTOR = 0.0207 


Figure 35. Decays of Blade Response with an Applied Root Loading of 2879N 
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LOSS FACTOR r] 


0.12 



Figure 36. Damping at Fan Blade Root 
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7-24 
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000 


AXIAL LOAD — N 


Figure 37. Variation of Root Vibratory Stress with Root Loading 
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RESPONSE FREQUENCY 




Figure 39a, b. Schematic of Shroud Rubbing 



nrm 

Figure 40. Analytical Model for Shroud Rubbing 
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FIXTURE/VACUUM CHAMBER 



Figure 41. Shroud Loading Device for Investigation of Fan Blade Shroud Damping Mechanism 



Figure 42. Shroud Loading Mechanism 
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EXCITATION SYSTEM RECORDING SYSTEM 



Figure 43. Block Diagram of Basic Excitation and Recording Systems Used to Obtain 
Strain and Acceleration Response Data from Test Items Excited by the 
Electro-dynamic Shaker 
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INSTRUMENTATION PORT 


ii SHROUD LOADING 
I DEVICE i 


SHROUD 

ACCELEROMETERS 


ALUMINUM EIXTURE 
(VACUUM CHAMBER) 


EVACUATION PORT 


LDED AT ROOT 


Figure 44. Test Set-Up for Shroud Damping Investigation 
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NODAL 



FREQUENCY, Hz 


FREQUENCY, Hz 

RESPONSE DESCRIPTION 

85 

NONLINEAR — MODIFIED FIRST FREE BLADE BENDING MODE 

140 

NONLINEAR — HIGH FREQUENCY RESPONSE 

170 

NONLINEAR — HIGH FREQUENCY RESPONSE 

260 

ABOVE SHROUD BLADE BENDING (FLAP) MODE (4 =0.003) 

300 

NONLINEAR — MODIFIED SECOND FREE BLADE BENDING MODE 

475 

ABOVE SHROUD BLADE SECOND BENDING MODE {V - 0.003) 

510 

NONLINEAR — HIGH FREQUENCY RESPONSE 

760 

NONLINEAR — HIGH FREQUENCY RESPONSE 

880 

ABOVE SHROUD BLADE TORSION (4 = 0.001) 


Figure 45. Fan Blade Response Characteristics from a 2g Level Exploratory Sweep — 
Shrouds Loaded with 102N 
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m 

T 

(0 

i 


n i 1 I i i I I 1 1 [ I I 

80 100 200 300 400 500 1000 12001400 2000 

FREQUENCY, Hz ^ 


FREQUENCY, Hz 

RESPONSE DESCRIPTION 

93 

NONLINEAR — MODIFIED FIRST FREE BLADE BENDING MODE 

278 

ABOVE SHROUD BLADE BENDING (FLAP) MODE 

1109 

ABOVE SHROUD BLADE TORSION 

1242 

NONLINEAR — MODIFIED SECOND BENDING MODE 

1442 

ABOVE SHROUD SECOND BENDING MODE 


Figure 46. Fan Blade Response Characteristics from a 2g Level Exploratory Sweep — 
Shrouds Loaded with 150 N 
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Figure 47. First Fiap Mode Response to a Siow Frequency Sweep at Constant Input 
Acceleration Levels 
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ASMT STRESS — MPa (pk) 


■ — r= 

ZZ NORMAL SHROUD LOAD: 150 N 


MODE: FIRST FLAP 



TIME (0.1 second/cm) 


Figure 48. Damping Due to Rubbing at Shroud for a Typical Fan Blade — 
Variation of Response Decay with Initial Stress 
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TIME (0.1 second/cm) 


Figure 49. Damping Due to Rubbing at Shroud for a Typicai Fan Biade — Variation of 
Response Decay with Shroud Load and initiai Stress 
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ASMT ROOT STRESS — MPa (pk) 



Figure 50. Shroud Damping Effects on Blade Response (First Flap Mode -280 Hz) 
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Figure 52. Wear Areas on Shroud Load Platens 
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Figure 53. Friction Test Assembly (Upper Arm Raised) 

81-78-B 81-3-7-50 
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FREQUENCY; 140 Hz 
INPUT ACCELERATION: 5g pk 


BOLT STRAIN (FORCE) — NO MOTION 


LOWER PIECE ACCELERATION 


250 Hz LP FILTERED UPPER PIECE 
ACCELERATION 


PIECE ACCELERATION 


STRAIN (FORCE) 


STRAIN — LOWER PIECE 
ACCELERATION (FILT) LOOP 


Figure 55a. Typical Signals from Acceleration and Force Transducers 



FREQUENCY: 140 Hz 
INPUT ACCEL: 5g pk 
HORIZ AXIS; 0.0127 mm/div 
VERT AXIS; 115.2 N/div 


271 N ) 

494 N [• APPLIED NORMAL LOAD 
672 N ) 


FREQUENCY; 280 Hz 
INPUT ACCEL; 20 g pk 
HORIZ AXIS: 0.032 mm/div 
VERT AXIS: 115.2 N/div 


Figure 55b. Typical Measured Friction Force — Slip Loops 
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FRICTION COEFFICIENT 


KEY; 


KEY; FREQ, Hz/VEL, mm/sec 


NORMAL LOAD 


o 

80/33 

O 

271 N 

• 

1 40/34 

□ 

494 N 

□ 

140/56 

A 

672 N 

■ 

280/56 



A 

280/112 



A 

420/74 





NORMAL LOAD — N 


FREQUENCY — Hz VELOCITY — mm/sec (pK) 


00 

T 

u 

I 
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Figure 56 Variation of Equivalent Friction Coefficient During Sinusoidal Motion 
With Normal Load, Frequency and Maximum Relative Velocity 




Figure 57a. Schematic of Turbine Blade and Platform Damper 



mm 


Figure 57b. Analytical Model for Platform Damping 
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Figure 59. Turbine Blade Test Piece Geometry 
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Figure 60. Damper Stress Response Signals for Varying Input Acceleration Levels 



Figure 61. Baseline Conditions for Turbine Blade 
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Figure 62. Platform Damping Test Results (Large Damper) 
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Figure 63. Platform Damping Test Results (Small Damper) 
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Figure 64. Turbine Blade Platform Damping Characteristics (Small Damper) 
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